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Abstract 

The  purpose  of  this  research  is  to  improve  the  ability  to  image  Special  Nuclear 
Material  (SNM)  using  a  High  Purity  Germanium  Double  Sided  Strip  Detector 
(HPGeDSSD)  for  use  in  the  field  of  nuclear  nonproliferation.  When  used  for  gamma 
spectrometry,  HPGe  detectors  can  efficiently  distinguish  various  isotopes  because  of  their 
good  energy  resolution.  The  detector  used  for  this  research  is  a  liquid  nitrogen  cooled 
germanium  detector  with  five  charge  collection  strips  on  each  face.  The  strips  are 
orthogonal  to  each  other,  creating  a  five  by  five  array  of  data  collection  pixels.  When  the 
25  mm  by  25  mm  detector  array  is  coupled  to  a  converging  collimator  an  imaging  system 
is  achieved  which  provides  the  ability  to  image  SNM  sources  which  are  larger  than  the 
detector  crystal.  Eor  this  work,  an  ideal  converging  collimator  was  designed,  built  and 
employed.  The  new  collimator  provided  proper  spatial  resolution  of  all  sources  tested 
while  still  providing  sufficient  efficiency  for  imaging.  Another  system  improvement  is 
the  addition  of  an  OR  gate  which  provides  proper  timing  information,  resulting  in  the 
ability  to  image  using  three  strip  events  for  the  first  time.  A  new  image  processing  code 
was  developed  which  employs  the  ability  to  use  three  strip  events  while  reducing  the  time 
required  to  produce  an  image  from  approximately  five  hours  down  to  less  than  a  minute. 
Additionally,  the  new  code  demonstrates  the  ability  to  distinguish  between  multiple 
isotope  sources  in  a  single  image  through  the  simultaneous  use  of  multiple  energy 
windows.  The  system  improvements  developed  during  this  work  combine  to  improve  the 
overall  efficiency  and  flexibility  of  the  imaging  system.  The  capability  to  distinguish  a 
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point  source  from  a  distributed  source,  with  proper  spatial  resolution,  was  demonstrated. 
This  ability  is  the  primary  goal  for  developing  this  teehnology. 
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IMPROVED  MULTINUCLIDE  IMAGING  OL  SPECIAL  NUCLEAR  MATERIAL 
USING  A  HIGH  PURITY  GERMANIUM  ORTHOGONAL  STRIP  DETECTOR 


I.  Introduction 


Motivation 

The  Defense  Threat  Reduction  Agency  (DTRA),  as  part  of  fulfilling  their 
mission  of  safeguarding  American  interests  from  weapons  of  mass  destruction,  has  been 
tasked  to  develop  quality  tools  and  services  for  detecting  and  monitoring  special  nuclear 
material  (SNM).  On  site  inspections  of  countries  that  are  signatories  of  various  arms 
control  treaties  are  a  major  part  of  DTRA’ s  SNM  monitoring  mission.  Some  examples 
are  the  inspections  and  monitoring  operations  at  facilities  in  Russia,  Ukraine,  Kazakhstan 
and  Belarus  in  support  of  the  Strategic  Arms  Reduction  Treaty  (START). 

The  purpose  of  these  inspections  is  to  verify  and  determine  the  presence  and 
quality  of  SNM.  Very  often  the  SNM  is  contained  within  nuclear  weapons  which  the 
countries  are  permitted  to  possess  in  accordance  with  arms  control  treaties.  In  these 
cases,  a  major  role  of  the  inspections  is  to  ensure  that  the  actual  weapon  SNM  has  not 
been  removed  and  replaced  with  a  spoof  material  designed  to  mimic  the  actual  weapon 
pit.  Eor  this  work,  the  SNM  is  assumed  to  be  weapons  grade  plutonium.  Proper 
verification  of  the  SNM  size,  shape  and  isotopic  composition  can  help  determine  the 
presence  of  a  possible  spoof  weapon. 

The  weapon  storage  configuration  for  inspection  is  assumed  to  involve  the 
weapon  remaining  in  a  storage  cask.  The  radius  of  the  cask  is  assumed  to  be  71  cm  [15]. 
Due  to  treaty  limitations,  only  passive  inspection  techniques  are  allowed.  This  limits  the 
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radiation  available  for  measurement  to  that  emitted  by  the  weapon  pit  with  enough 
energy  to  penetrate  the  eask.  Owing  to  attenuation,  only  gamma  radiation  above  200  keV 
is  energetie  enough  to  have  a  measurable  flux  that  ean  be  deteeted  on  the  outside  of  the 
eask  [14],  The  information  required  to  verify  the  possibility  of  a  spoof  weapon  ineludes 
the  gamma  speetrum  of  the  weapons  grade  plutonium  (WGPu)  along  with  the  basie  size 
and  shape  of  the  pit.  Deteetors  eurrently  used  for  on  site  inspeetions  are  not  eapable  of 
distinguishing  whether  a  weapon  eontains  SNM  in  the  form  of  a  point  souree  or  a 
distributed  souree.  They  also  eannot  distinguish  WGPu  from  plutonium  having  less 
fissile  Pu-239  eontent. 

Background 

Previous  researeh  at  the  Air  Foree  Institute  of  Teehnology  demonstrated 
the  eapability  of  differentiating  WGPu  from  reaetor  grade  plutonium  using  the  same 
HPGeDSSD  used  for  this  work  [14].  In  addition,  the  use  of  a  eollimator  made  it  possible 
to  distinguish  a  point  souree  from  a  distributed  souree  loeated  at  the  eollimator  front  faee. 
The  previous  eollimator  utilized  a  set  of  parallel  holes.  This  limited  the  size  of  the  souree 
to  that  of  the  deteetor  erystal,  25  mm  by  25  mm.  Further  analysis  of  the  previous  design 
also  showed  that  the  field  of  view  for  eaeh  eollimator  hole  on  the  deteetor  was  larger  than 
a  single  pixel.  This  led  to  an  inaoeurate  spatial  resolution  of  the  souree.  Using  a 
eonverging  eollimator,  one  should  be  able  to  determine  the  basie  shape  and  size  of  a 
larger  souree.  Proper  spatial  resolution  ean  be  maintained  by  ensuring  that  eaeh 
eollimator  hole  field  of  view  eovers  only  a  single  pixel. 
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The  detector  is  smaller  than  conventional  HPGe  detectors,  which  results  in  lower 
detection  efficiencies  for  higher  gamma  energies.  However,  as  a  significant  trade-off,  the 
orthogonal  strip  array  provides  the  position  sensitivity  required  to  image  SNM  contained 
within  a  weapon  pit.  The  physical  dimensions  of  a  weapon  pit  are  discussed  in  Chapter 

2. 

Strip  Detector  Imaging;  Position  sensitivity  is  a  crucial  requirement  for  the 
development  of  a  detector  system  capable  of  producing  an  image.  The  HPGe  double 
sided  strip  detector  design  provides  the  position  sensitivity  required  for  imaging  without 
sacrificing  the  energy  resolution  required  for  differentiating  WGPu.  The  previous  work 
conducted  with  this  detector  proved  that  a  collimated  system  can  produce  the  data 
required  to  develop  a  source  image.  The  previous  source  size  limitation  can  be  overcome 
by  mating  the  strip  detector  to  a  converging  hole  collimator.  A  collimator  is  used  to 
reduce  incident  photons  to  a  known  range  of  incident  angles.  While  the  absolute 
efficiency  of  the  detector  is  greatly  reduced  by  the  reduction  of  incident  photons  on  the 
detector,  it  provides  the  ability  to  image  a  source  larger  than  the  detector  strip  array. 
Figure  1  shows  the  coupling  of  a  detector  with  a  converging  collimator  to  create  an 
imaging  system. 
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Isotopic  Source  Absorbing  Converging  Collimator  Detector 

Figure  1:  Converging  collimator  imaging  system 

Imaging  Device:  The  deviee  used  for  imaging  employs  an  HPGe  orthogonal 
strip  detector  with  25  pixels  in  a  five  by  five  array,  as  shown  in  Figure  2.  Each  pixel 
measures  5mm  by  5mm.  Due  to  hardware  limitations,  only  a  four  by  four  array  is 
utilized  for  imaging  during  this  work.  However,  the  collimator  is  designed  for  use  with 
the  entire  five  by  five  array  to  allow  the  flexibility  of  choosing  which  four  by  four 
combination  is  used. 
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In  order  to  produce  an  image  with  proper  spatial  orientation,  a  collimator  is 
required  to  ensure  that  the  counted  photons  originated  from  a  particular  source  volume. 
The  imaging  collimator  for  this  work  is  60  mm  by  60  mm  by  190  mm  and  is  made  of 
CerroBEND,  a  lead  bismuth  alloy.  It  has  a  5  by  5  array  of  25  square  holes,  each  2.9  mm 
wide.  Although  the  holes  have  a  constant  dimension  along  their  entire  length,  they  are 
arranged  in  a  converging  pattern  to  allow  imaging  of  a  source  larger  than  the  detector 
crystal.  Proper  alignment  of  the  collimator  with  the  detector  and  source  ensures  that  each 
collimator  hole  produces  a  source  view  in  only  one  detector  pixel,  thus  providing  proper 
spatial  resolution  of  the  source.  Figure  3  shows  the  collimator  fields  of  view 
superimposed  on  the  detector  crystal  as  calculated  using  the  collimator  design  code 
presented  in  Appendix  A.  The  black  holes  represent  the  configuration  on  the  rear  of  the 
crystal.  The  grey  holes  represent  the  configuration  on  the  front  of  the  crystal. 


Figure  3:  Collimator  fields  of  view  on  the  detector  faces 

Data  Acquisition:  The  output  signals  from  the  HPGeDSSD  preamplifiers  are  fed 
to  a  series  of  Digital  Gamma  Finder  (DGF)  modules.  The  modules  process  the  digital 
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pulses  and  produce  an  output  fde  consisting  of  module  identification,  the  charge 
collection  strip  hit  patterns,  timestamps  and  channel  numbers  for  each  event  acquired.  In 
the  previous  work  using  the  HPGeDSSD  system,  it  was  not  possible  to  stop  and  start  all 
of  the  DGFs  at  the  same  time.  This  made  the  assignment  of  events  coincident  on 
opposite  faces  of  the  detector  difficult.  The  difficulty  was  in  the  inability  to  obtain  proper 
timing  information  because  once  a  module  filled  its  first  level  of  memory,  referred  to  as  a 
spill,  it  would  download  its  data  to  the  output  file.  Meanwhile,  the  other  modules  would 
continue  acquiring  data  until  their  first  level  memory  was  full,  at  which  time  they  would 
download  their  data  and  begin  a  new  spill.  Since  the  timestamp  recorded  for  an  event  is 
counted  from  the  time  the  current  spill  began,  after  the  first  spill  of  information  was 
downloaded  each  module  would  be  recording  different  times  for  events  which  were 
actually  coincident  in  time.  Therefore,  the  inability  to  start  and  stop  all  modules  at  the 
same  time  degraded  the  efficiency  of  the  system.  This  problem  should  be  overcome  by 
the  addition  of  an  input  logic  module  that  will  maintain  consistent  timing  information 
between  all  modules.  This  should  result  in  the  ability  to  start  and  stop  all  module  data 
acquisitions  at  the  same  time. 

Objective 

The  focus  of  this  thesis  is  to  develop  an  improved  imaging  system  using  an  HPGe 
double-sided  strip  detector  (HPGeDSSD).  The  first  improvement  is  the  design  and 
construction  of  a  converging  hole  collimator  which  maintains  proper  spatial  resolution  of 
the  source  without  sacrificing  the  efficiency  of  the  system.  This  includes  being  able  to 
image  a  source  similar  in  size,  shape  and  configuration  to  that  encountered  during  on  site 
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weapon  inspections.  The  convergence  angles  of  the  collimator  holes  allow  imaging  of  a 
source  consistent  in  size  and  shape  with  that  encountered  during  on  site  inspections. 

Using  a  technique  similar  to  a  simple  Anger  camera  [1],  the  25  detector  pixels  are 
capable  of  producing  an  image  showing  the  basic  shape  and  size  of  a  distributed  source. 
This  provides  the  ability  to  differentiate  a  distributed  source  from  a  point  source. 

An  additional  system  improvement  includes  a  more  efficient  image  processing 
code.  The  previous  imaging  code  required  approximately  five  hours  of  data  processing  to 
produce  a  single  image.  Shortening  this  to  a  few  minutes  would  make  the  system  useful 
for  on  site  weapons  inspections.  The  previous  code  was  also  limited  to  imaging  across 
only  a  single  energy  window.  A  code  which  can  produce  images  using  multiple  energy 
windows  would  provide  greater  flexibility  for  data  acquisition  and  analysis. 

The  final  system  improvement  involves  the  addition  of  an  input  logic  module 
which  should  force  all  DGF  modules  to  start  and  stop  at  the  same  time.  This  should 
result  in  proper  relative  timing  of  all  acquired  pulses  and  thus  allow  proper  assignment  of 
coincident  pulses. 

The  photon  energies  used  for  imaging  must  be  high  enough  to  escape  a  weapon 
cask  with  a  measurable  flux.  Since  the  efficiency  of  a  germanium  detector  decreases 
with  increasing  photon  energy,  the  energy  of  the  imaging  photons  must  be  in  a  range 
which  balances  acquisition  with  detector  efficiency.  The  imaging  of  SNM  in  a  weapon 
pit  can  be  done  by  observing  the  gamma  radiation  from  the  dominate  isotope,  Pu-239, 
which  has  prominent  peaks  at  413.71  keV  and  375.05  keV.  Both  of  these  peaks  should 
be  high  enough  in  energy  to  escape  from  the  system  and  still  be  low  enough  in  energy  to 
be  efficiently  detected  and  imaged.  The  413  keV  peak  has  a  yield  of  0.00147%  and  the 
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375  keV  peak  has  a  yield  of  0.00155  %.  In  an  operational  system  one  or  both  of  these 
peaks  are  proposed  to  be  used  for  imaging  the  pit.  For  mueh  of  this  work,  the  SNM  will 
be  simulated  using  a  neutron  activated  Sn  source.  The  activated  Sn  has  a  prominent  peak 
at  391.73  keV  and  427.89  keV. 


Scope 

This  thesis  work  began  with  initial  analysis  of  the  collimator  design  using  the 
Mathematica  code  in  Appendix  A.  Collimator  Design  Code  Using  the  results  of  the 
design  code,  several  collimators  were  made  to  determine  detector  response  to  varying 
hole  sizes.  This  also  led  to  a  test  to  determine  the  feasibility  of  using  holes  located  at  the 
edge  of  the  detector  pixels  in  order  to  increase  the  spatial  resolution  of  the  system.  These 
tests  were  also  used  to  characterize  the  response  of  gamma  radiation  incident  at  several 
locations  throughout  a  pixel,  as  well  as  determine  the  settings  to  ensure  that  each  strip 
had  the  proper  energy  calibration.  The  results  of  the  collimator  tests  led  to  the  design  of  a 
25  hole  converging  collimator.  It  was  then  shown  that  mating  the  HPGeDSSD  to  the 
converging  hole  collimator  produces  the  ability  to  image  both  point  and  distributed 
gamma  ray  sources.  The  raw  data  was  processed  using  the  MATLAB  code  shown  in 
Appendix  B.  Image  Processing  Code,  resulting  in  a  four  by  four  pixel  image  of  the 
source,  a  surface  plot  of  the  acquired  image  and  a  histogram  showing  the  number  of 
counts  acquired  at  energies  within  user  specified  energy  windows.  The  user  can  choose 
to  use  one  or  two  energy  windows. 
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Assumptions 


In  order  to  best  simulate  an  operational  environment,  several  assumptions  have  to 
be  made.  The  first  is  an  assumption  that  the  SNM  to  be  imaged  is  symmetrieal,  sinee  the 
image  is  only  two  dimensional.  However,  while  this  assumption  must  be  made  during 
inspections,  this  project  uses  a  source  which  is  not  always  symmetrical  for  demonstrating 
the  ability  of  the  HPGeDSSD  to  image.  The  second  assumption  is  that  the  maximum 
amount  of  time  allowed  to  image  weapon  pits  during  inspections  will  be  30  minutes, 
however  this  work  does  not  impose  restrictions  on  the  counting  times.  The  third 
assumption  is  that  the  SNM  located  within  the  pit  is  plutonium.  Finally  it  is  assumed  that 
during  the  inspection  process,  the  SNM  cannot  be  removed  from  the  71  cm  radius  steel 
cask,  which  requires  that  only  passive  inspections  be  utilized. 


General  Approach 

Calibration  of  each  detector  strip  to  the  same  energy  scale  is  necessary  to  produce 
an  image  of  a  specific  isotope  contained  within  a  radioactive  source.  Proper  alignment  of 
the  converging  collimator  must  then  be  verified  in  order  to  produce  the  desired  image. 
Proper  settings  for  the  data  acquisition  components  help  to  produce  raw  data  which  is 
processed  through  a  MATLAB  code  to  produce  an  image. 


Sequence  of  Presentation 

Chapter  2  gives  the  reader  an  introduction  to  the  physics  of  germanium  detectors, 
the  pulse  processing  required,  an  introduction  to  strip  detectors,  and  finally  how  the  strip 
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detector  is  used  to  create  an  image.  Chapter  3  describes  the  detector,  processing 
electronics,  controlling  software,  collimator,  image  processing  code,  and  radiation 
sources  used  during  the  experiment.  Chapter  4  details  the  experimental  and 
computational  procedures  employed  during  the  project.  Chapter  5  presents  the  results 
obtained  during  the  project,  and  Chapter  6  finishes  with  conclusions  and 
recommendations  for  further  work. 
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II,  Theory 


Chapter  Overview 

This  chapter  introduces  the  prineiples  of  germanium  deteetors,  pulse 
proeessing,  and  strip  deteetors.  These  prineiples  are  the  foundation  for  the  experimental 
setup  and  analysis,  whieh  are  discussed  in  Chapters  4  and  6  respeetively. 


Germanium  Detectors  [6] 

Germanium  deteetors,  like  all  semieonductor  deteetors,  operate  by  the  eolleetion 
of  eleetron-hole  pair  charges  produeed  by  ionizing  radiation  in  the  deteetor  material.  Of 
the  various  ways  gamma  rays  can  produce  electron-holes  pairs  in  the  germanium  erystal, 
only  two  interaetion  meehanisms  have  any  signifieanee  for  this  work:  photoeleetrie 
absorption  and  Compton  seattering.  Both  proeesses  result  in  the  production  of  electron- 
hole  pairs  within  the  germanium.  Under  the  influenee  of  an  applied  eleetrie  field,  the 
negatively-eharged  eleetrons  will  drift  through  the  erystal  following  the  lines  of  foree 
produeed  by  the  eleetrie  field.  The  hole,  representing  a  net  positive  eharge,  will  also  tend 
to  move  in  the  eleetrie  field,  but  in  a  direetion  opposite  that  of  the  eleetron.  The  direetion 
of  motion  for  eaeh  is  approximately  perpendicular  to  the  detector  electrical  contacts.  The 
motion  of  both  of  these  eharges  eontributes  to  the  observed  eonductivity  of  the 
germanium,  resulting  in  a  eharge  pulse  being  colleeted  on  the  deteetor  eleetrical  eontaets. 
The  amplitude  of  the  resulting  eharge  pulse  is  proportional  to  the  amount  of  energy 
deposited  by  the  gamma  photon. 
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The  number  of  electron-hole  pairs  produced  by  a  gamma  photon  is  a  function  of 
the  stopping  power  (dE/dX),  and  the  energy  required  to  produce  an  electron-hole  pair. 

The  number  of  electron-hole  pairs  that  contribute  to  the  output  signal  is  determined  by 
the  charge  collection  efficiency,  which  is  dependent  on  the  carrier  mobility,  lifetime  of 
the  carriers,  and  the  collection  time.  The  collection  time  decreases  with  an  increasing 
field  gradient  within  the  detector  volume. 

Efficient  charge  collection  is  achieved  by  ensuring  the  applied  bias  voltage 
produces  full  charge  carrier  depletion  from  the  entire  active  volume  of  the  detector. 
Achieving  full  depletion  reduces  the  detector  capacitance  and  maximizes  the  signal-to- 
noise  ratio,  making  it  a  preferred  configuration  in  most  applications.  The  depletion  depth 
becomes  very  important  for  application  at  gamma  ray  energies  from  200  to  700  keV.  The 
greater  the  depletion  region,  the  better  the  detector  will  be  at  stopping  higher  energy 
gammas.  The  thickness  of  the  depletion  region  d  for  a  planar  detector  configuration  is 
given  by  [6] 


"levy 


(1) 


where; 

V  is  the  reverse  bias  voltage, 

N  is  the  net  impurity  concentration  in  the  bulk  semiconductor  material, 
e  is  the  electronic  charge,  and 
s  is  the  dielectric  constant  (16  So  for  germanium). 
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Therefore  at  a  fixed  bias  voltage,  the  only  way  to  inerease  the  depletion  depth  is 
to  lower  the  impurity  concentration.  For  the  detector  used  in  this  research,  this  was 
achieved  by  using  high  purity  germanium,  where  the  impurity  concentrations  are 
approximately  10^°  atoms/cm^.  This  allows  for  a  depletion  region  of  about  13mm  at 
lOOOV,  thus  fully  depleting  the  9  mm  thick  crystal  used  for  this  work. 

Figure  4  shows  the  configuration  of  a  planar  HPGe  detector  fabricated  out  of  p- 
type  germanium  [6]. 


Ge  disk 


Figure  4:  Planar  HPGe  detector  configuration  [6] 


Planar  HPGe  detectors  have  electrical  contacts  on  two  flat  surfaces.  The  n+  contact  is 
formed  by  lithium  evaporation,  and  the  p+  contact  is  formed  by  ion  implantation  of  boron 
atoms.  When  the  junction  is  reverse  biased  a  depleted  region  develops.  Typically,  HPGe 
detectors  are  operated  as  fully  depleted  detectors.  Reverse  biasing  of  the  detector 
required  that  a  positive  voltage  be  applied  to  the  n+  contact  of  the  detector  and  the  p+ 
contact  be  maintained  at  a  relative  ground.  When  planar  germanium  detectors  are  fully 
depleted  and  are  operated  with  a  large  over  voltage,  the  electric  field  can  be  assumed  to 
be  uniform  from  one  contact  to  the  other.  This  allows  the  electrons  and  holes  to  be 
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treated  as  drifting  under  the  influence  of  a  constant  electric  field  through  the  entire 
volume  of  the  detector. 

The  most  important  parameters  for  any  gamma  spectrometer  are  the  efficiency  and 
energy  resolution.  The  efficiency  determines  the  number  of  counts  that  can  be  collected 
over  a  given  time,  while  the  energy  resolution  is  a  measure  of  the  ability  to  distinguish 
gamma-rays  with  similar  energies.  For  germanium  detectors  the  resolution  is  expressed 
as  the  full  width  at  half  maximum  (FWHM)  for  a  particular  full  energy  peak. 

Pulse  shaping/processing  [17] 

Collection  of  the  electron-hole  pairs  produces  a  charge  pulse.  Each  charge  pulse 
is  processed  by  a  preamplifier  which  is  built  into  the  detector  cryostat  housing.  The 
magnitude  of  the  voltage  pulses  sent  from  the  preamplifier  to  the  pulse  processing 
electronics,  is  dependent  on  the  amount  of  charge  collected  and  therefore  dependent  on 
the  energy  that  was  deposited  and  collected. 

When  the  pulse  is  sent  from  the  preamplifier  to  the  pulse  processing  electronics, 
the  pulse’s  waveform  is  digitized  in  an  analog  to  digital  converter  (ADC).  Next  the 
signal  is  passed  to  a  real  time  processing  unit  (RTPU)  where  a  digital  filter  is  applied  to 
shape  and  amplify  the  pulse.  The  filter  used  is  a  trapezoidal  filter;  this  filter  typically  has 
a  fiat  top,  which  covers  the  rise  time  of  the  pulse  that  makes  the  pulse  height 
measurement  less  sensitive  to  the  variations  of  the  signal  shape.  An  example  of  a 
trapezoidal  fdter  applied  to  an  ADC  output  pulse  is  shown  in  Figure  5. 
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Figure  5;  Trapezoidal  filter  applied  to  a  pre-amplifier  pulse  [17] 

The  rise  time  L,  and  the  flat  top  G  of  the  filter  ean  be  controlled  by  the  user 
through  the  XIA  controlling  software  to  maximize  the  energy  resolution. 


Strip  Detectors 

A  germanium  orthogonal  strip  detector  combines  the  energy  resolution  of  a  high 
purity  germanium  detector  with  the  ability  to  produce  an  image.  Historically  the  high 
energy  resolution  and  imaging  capabilities  of  the  germanium  strip  detectors  have  been 
exploited  through  astrophysics  applications  [8].  Germanium  strip  detectors  are  similar  to 
planar  HPGe  detectors,  except  that  single  contacts  to  collect  the  charges  on  each  side  are 
replaced  with  discrete  strips  of  contacts  which  are  orthogonal  to  each  other.  Figure  6 
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shows  the  configuration  of  the  charge  collection  strips  for  a  five  by  five  strip  detector, 
like  that  used  for  this  work. 


+Hf 


Figure  6:  Electrode  strip  configuration  of  a  strip  detector  [8] 

With  a  standard  FlPGe  planar  detector  the  charge  is  typically  read  out  from  one  of 
the  two  contacts,  usually  the  n+  contact  which  collects  the  electrons.  However,  with  a 
germanium  strip  detector  all  of  the  contacts  are  read  out  to  individual  preamplifiers. 

Since  each  interaction  within  the  germanium  must  produce  a  hole  and  an  electron,  at  least 
one  of  the  strips  on  each  side  must  collect  the  charges.  This  ability  to  collect  electrons 
and  holes  on  opposite  sides  of  the  germanium  is  what  provides  the  means  for  producing 
an  image. 
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Imaging  Using  a  Strip  Detector 


The  intersections  of  the  orthogonal  electrode  strips  which  collect  the  electrons  and 
holes  produce  a  five  by  five  array  of  pixels.  The  ability  to  distinguish  which  strips  collect 
the  charge  from  a  particular  ionizing  interaction  provides  a  way  to  obtain  two- 
dimensional  position  localization  of  the  interaction. 

Following  the  recommendation  from  previous  imaging  work  performed  using  this 
system,  I  proposed  using  two  categories  of  image  pixels.  First  are  the  “center  pixels” 
which  correspond  to  collimator  holes  centered  on  one  of  the  25  detector  pixels  formed  by 
the  intersection  of  a  front  and  back  strip.  Second  are  the  “edge  pixels”  which  correspond 
to  collimator  holes  centered  on  the  intersection  of  two  strips  on  one  face  and  centered  on 
a  strip  on  the  opposite  face.  Figure  7  shows  the  relative  locations  of  the  two  categories  of 
pixels  superimposed  onto  the  25  detector  pixels.  In  theory,  the  use  of  edge  pixels  would 
increase  the  spatial  resolution  of  the  system  by  narrowing  the  origin  of  the  gamma  ray 
down  to  a  smaller  area. 


^  -  Center  pixels 
Q  -  Edge  pixels 


Figure  7:  Center  and  edge  pixel  positions 
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In  order  for  an  event  on  a  eenter  pixel  to  be  used  in  the  formation  of  the  image,  a 
front  strip  and  a  baek  strip  must  both  reeord  an  event  that  oeeurred  at  the  same  energy 
and  at  the  same  time.  To  be  eonsidered  the  same  energy  two  pulses  must  be  within  one 
FWHM  of  a  full  energy  event  peak.  To  be  eonsidered  the  same  time,  the  two  pulses  must 
fall  within  some  eoineidenee  window  of  time.  For  this  projeet  the  time  eoineidenee 
window  was  set  to  100  ns.  This  was  optimized  experimentally. 

Therefore,  if  a  front  strip  and  a  baek  strip  have  a  eoineident  event  in  energy  and 
time,  the  loeation  of  the  interaetion  ean  be  inferred  as  the  plaee  where  the  two  strips 
interseet.  This  is  referred  to  as  a  two  strip  event.  These  events  are  the  result  of  a  gamma 
depositing  energy  in  only  one  pixel  of  the  deteetor.  For  example,  if  a  eharge  eolleetion 
strip  on  the  front,  F2,  and  a  eharge  eolleetion  strip  on  the  baek,  R3,  both  reeord  events 
eoineident  in  time  and  energy,  then  the  interseetion  of  the  two  strips  is  the  loeation  of  the 
event,  as  shown  in  Figure  8. 


• 

F2 

R3 

Figure  8:  Position  loealization  for  a  eenter  pixel 
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In  order  for  a  three  strip  event  on  an  edge  pixel  to  be  used  in  the  formation  of  the 
image,  one  strip  on  either  the  front  or  baek  must  reeord  a  full  energy  peak  for  the  energy 
being  sought.  In  addition,  the  energy  peaks  on  two  adjaeent  strips  on  the  opposite  side  of 
the  deteetor  must  sum  up  to  an  energy  that  is  within  one  FWHM  of  the  full  energy  event 
peak  of  the  original  strip.  All  three  peaks  must  fall  within  the  time  eoineidenee  window, 
just  as  with  the  eenter  pixel  two  strip  events.  For  example,  suppose  a  eharge  eolleetion 
strip  on  the  baek,  R3,  reeords  an  event.  In  order  to  eount  as  an  edge  pixel  event,  there 
must  be  two  adjaeent  strips  on  the  front,  say  F3  and  F4,  whose  eombined  event  is 
eoineident  in  both  time  and  energy  with  the  event  reeorded  on  R3.  This  would  result  in 
the  assignment  of  the  event  to  the  edge  pixel  shown  in  Figure  9. 


Figure  9;  Position  loealization  for  an  edge  pixel 

In  order  for  sueh  an  event  to  be  reeorded,  one  of  two  things  would  have  to  oeeur. 
First,  a  gamma  eould  undergo  a  photoeleetrie  absorption  very  elose  to  the  gap  between 
the  adjaeent  strips.  If  elose  enough  to  the  gap,  some  eharge  would  be  eolleeted  by  one 
strip  and  the  rest  by  the  adjaeent  strip.  The  second  scenario  involves  single  or  multiple 
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Compton  scatters  which  send  the  reduced  energy  gamma  into  the  adjacent  pixel  where  it 
would  undergo  either  a  photoelectric  absorption  or  more  Compton  scatters  followed  by 
photoelectric  absorption.  In  both  scenarios  the  adjacent  strips  would  show  a  charge 
sharing  of  the  full  energy  event  recorded  on  the  opposite  side  of  the  crystal. 

As  discussed  in  later  chapters,  the  use  of  collimator  holes  producing  edge  pixels 
proved  to  be  impracticable  for  gamma  imaging.  However,  if  an  event  shares  charge 
between  two  adjacent  strips  and  one  opposite  face  strip,  as  discussed  for  the  edge  pixels, 
it  can  still  be  used  for  imaging  if  proper  timing  information  is  acquired.  If  the  charge 
collection  on  adjacent  strips  occurs  within  the  time  resolution  of  the  system,  namely  25 
ns,  the  interactions  appear  to  occur  simultaneously  in  two  pixels  and  therefore  cannot  be 
used.  However,  if  the  timestamps  for  the  two  adjacent  strip  charge  collections  is 
separated  by  more  than  25  ns  the  interaction  timestamps  will  be  separated  by  at  least  one 
time  unit.  This  separation  in  time  allows  the  event  to  be  assigned  to  the  pixel  which 
recorded  the  interaction  first.  This  assignment  of  three  strip  events  to  a  single  pixel  is 
based  on  the  assumption  that  the  gamma  interaction  resulted  in  Compton  scattering  of  the 
gamma  into  the  adjacent  pixel.  The  time  required  for  the  gamma  to  scatter  into  the 
adjacent  pixel,  interact,  and  produce  a  signal  produces  the  time  difference  between  the 
interactions  recorded  in  the  two  pixels.  This  assumption  requires  that  the  incident  gamma 
be  well  collimated  to  ensure  that  the  first  Compton  scatter  interaction  occurs  in  the  pixel 
where  the  gamma  first  entered  the  detector. 

If  the  source  irradiating  the  detector  is  well  collimated,  and  the  location  of  the 
interaction  within  the  detector  is  known,  the  location  from  which  the  photon  emanated 


20 


can  also  be  known.  The  uneertainty  in  the  photon’s  loeation  of  origin  is  limited  by  the 
size  of  the  eollimator  field  of  view  on  the  souree. 


Image  Analysis  [13] 

The  image  analysis  for  this  projeet  involved  images  eonsisting  of  a  four  by  four 
array  of  grayseale  pixels.  Eaeh  pixel  eorresponds  to  one  pixel  formed  by  two  of  the 
orthogonal  deteetor  charge  eolleetion  strips.  Eaeh  event  recorded  by  the  deteetor  is 
analyzed  to  see  if  it  meets  the  speeified  requirements  to  be  eounted  as  an  image  event.  If 
the  image  requirements  are  met,  then  the  event  is  assigned  to  its  appropriate  pixel 
loeation.  When  all  events  have  been  proeessed,  the  resulting  four  by  four  matrix  is  used 
to  eonstruct  the  grayseale  image.  The  range  of  values  for  each  pixel  will  be  from  zero  to 
the  number  of  events  reeorded  for  that  partieular  pixel,  whieh  may  be  as  high  as  several 
thousand  depending  upon  the  source.  The  maximum  pixel  value  in  the  image  is  then 
assigned  a  value  representing  pure  white  (255).  A  eount  of  zero  eorresponds  to  blaek  (0), 
all  other  pixels  are  then  assigned  a  grayseale  value  based  upon  their  value  relative  to  the 
maximum  value  in  the  image. 

The  possibility  of  having  values  ranging  up  to  several  thousand  eounts 
eomplieates  the  imaging  proeess  beeause  there  are  only  256  distinguishable  brightness 
values  for  8-bit  grayscale  imaging.  These  values  eorrespond  to  an  8-bit  number  ranging 
from  0  to  255,  0  being  blaek  and  255  being  white.  Therefore,  if  the  image  eontains 
values  above  255  then  all  pixel  values  must  be  sealed  to  fit  within  the  0  to  255  range. 
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The  pixel  value  sealing  is  further  eomplieated  by  the  faet  that  only  about  20-30 
gray  levels  ean  be  visually  distinguished  by  the  human  eye.  However,  a  good  quality 
printing  system  ean  produee  more  than  30  gray  levels.  This  ean  lead  to  difficulty  in 
properly  distinguishing  between  different  pixel  values  in  an  image.  Providing  a  grayscale 
legend  for  an  image  can  greatly  aid  in  distinguishing  the  true  relative  brightness  between 
pixels.  However,  the  best  grayscale  image  analysis  is  performed  by  referring  to  the 
actual  pixel  values  rather  than  the  resulting  grayscale  image.  This  type  of  analysis  would 
be  quite  difficult  for  an  image  containing  several  thousand  pixels  but  it  works  quite  well 
for  the  16  pixel  images  produced  for  this  work. 

Weapon  Assumptions 

The  material  and  physical  dimensions  of  a  nuclear  weapon  are  key  components 
that  will  affect  the  ability  of  the  HPGeDSSD  to  produce  an  image  of  the  SNM.  The  SNM 
considered  for  this  project  is  assumed  to  be  Pu-239.  The  high  energy  portion  of  a  Pu-239 
spectrum  collected  with  a  germanium  detector  by  the  Idaho  National  Engineering  & 
Environmental  Eaboratory  is  shown  in  Eigure  10.  The  peaks  of  interest  for  this  project 
are  the  375  keV  (not  shown)  and  413  keV  peaks.  These  peaks  are  simulated  using  the 
391  keV  and  427  keV  peaks  from  a  neutron-activated  tin  source.  Eurther  discussion  of 
the  neutron-activated  tin  spectrum  can  be  found  in  Chapter  3. 
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Figure  10;  Pu-239  gamma  spectrum  from  an  HPGe  detector  [4] 


The  time  required  to  produce  an  image  depends  primarily  on  the  radiation  flux  at 
the  detector  that  are  to  be  imaged.  In  an  imaging  application  of  nuclear  weapons,  the  flux 
is  low  because  of  the  attenuating  materials  surrounding  the  SNM,  as  well  as  the  low 
activity  of  plutonium  due  to  its  long  half  life.  An  understanding  of  the  physical 
dimensions  of  the  pit  and  the  shielding  around  a  weapon  can  be  obtained  by  studying  a 
model  published  by  Steve  Fetter  in  “Detecting  Nuclear  Materials”  [3].  In  this 
publication,  Fetter  presents  what  he  refers  to  as  the  general  characteristics  of  an 
implosion- type  fission  weapon.  Figure  1 1  is  a  description  of  a  hypothetical  weapon 
presented  in  reference  [3]. 
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Fissile  core:  4  kg 
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High  explosive,  10  cm  thick 
Aluminum  case,  1  cm  thick 


Figure  1 1 :  Fetter  Model  of  a  hypothetical  implosion  weapon 
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III.  Equipment 


Chapter  Overview 

The  components  of  the  imaging  system  used  for  this  work  include  the  germanium 
strip  detector,  the  converging  hole  collimator,  two  digital  gamma  finder  (DGF)  data 
acquisition  cards,  the  input  logic  unit,  the  power  distribution  module  (PDM),  the 
CAMAC  crate,  the  crate  controller,  the  software  to  control  the  DGFs,  and  the  software  to 
produce  the  image.  A  schematic  of  the  entire  system  is  shown  in  Figure  12.  A  photo  of 
the  HPGeDSSD  imaging  system  electronics  configuration  is  shown  in  Figure  13. 


Figure  12:  HPGeDSSD  imaging  system  schematic 
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Figure  13;  Photo  of  HPGeDSSD  imaging  system  electronics 


Detector  [10]  [14] 

The  detector  used  for  this  project  was  an  ORTEC  Model  HPGeDSSD, 
Serial  number  42-WPAFB-Ol.  The  detector  is  an  HPGeDSSD  that  operates  at 
approximately  77K  using  liquid  nitrogen.  A  photo  of  the  detector,  dewar  and  collimator 
stage  is  shown  in  Figure  14. 
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Figure  14;  Detector,  dewar  and  collimator  stage 


The  detector  has  five  charge  collecting  strips  on  each  25  mm  by  25  mm  side.  It 
was  designed  and  specified  by  Mr.  Matt  Lange,  a  previous  AFIT  student.  The  strips  on 
opposing  sides  are  orthogonal  to  each  other,  which  allows  for  the  determination  of 
location  when  a  gamma  interaction  occurs  within  the  germanium  crystal.  This  effectively 
creates  25  pixels  on  the  detector.  Each  of  the  charge  collection  strips  is  25  mm  long  and 
5  mm  wide,  so  that  the  total  active  surface  area  of  the  detector  is  25  mm  by  25  mm.  The 
detector  has  an  active  depth  of  9  mm.  The  dimensions  of  the  detector  and  the  charge 
collection  strips  are  shown  in  Figure  15. 
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5mm 

Figure  15:  ORTEC  germanium  strip  deteetor  dimensions 


The  deteetor  is  seeured  inside  a  eylindrieal  eryostat  housing  whieh  allows  it  to  be 
eooled  to  77K.  The  housing  also  eontains  a  total  of  10  preamplifiers,  one  for  eaeh  eharge 
eolleetion  strip.  In  order  to  minimize  the  influenee  of  the  housing  on  the  ineident  gamma 
rays,  there  is  a  0.5  mm  thiek  beryllium  entranee  window,  as  well  as  0.3  pm  of  inaetive 
germanium  along  the  eenter  axis  of  the  housing.  The  deteetor  is  oriented  so  that  the  25 
mm  by  25  mm  faee  is  perpendieular  to  the  long  axis  of  the  eollimator  eenter  hole. 

In  order  to  optimize  the  eolleetion  of  eleetron-hole  pairs,  the  deteetor  is  operated 
at  a  +1000  V  bias.  The  eharge  eolleeting  strips  on  the  baek  faee  of  the  deteetor  are 
lithium  and  are  at  +1000  V  bias.  The  strips  on  the  front  faee  are  boron  and  are  at  about 
+1  V.  Eaeh  of  the  strips  is  eleetrieally  separated  from  one  another.  Sinee  there  is  a  1000 
V  potential  between  the  front  and  rear  strips,  the  two  sides  ean  be  treated  as  parallel 
plates,  whieh  foree  any  eleetrons  or  holes  to  the  eharge  eolleetion  strips. 
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DGF  Hardware  [17] 


For  this  project,  two  X-Ray  Instrumentation  Associates  (XIA)  DGF  model  4C, 
Revision  E-1  CAMAC  modules  were  used.  Their  serial  numbers  were  1240  and  1244. 
The  DGF-4C  is  a  four  channel  all-digital  waveform  acquisition  and  spectrometer  card. 
The  incoming  signals  from  the  detector  are  digitized  by  analog-to-digital  converters 
(ADC).  Digital  waveforms  are  analyzed  by  the  digital  signal  processor  (DSP)  for  pulse 
shape  analysis.  Then  the  waveforms,  timestamps,  and  the  pulse  shape  analysis  results  are 
read  by  the  host  system.  The  modules  can  process  up  to  200,000  counts  per  second  for 
all  four  channels  combined.  The  individual  modules  support  coincidence  spectroscopy 
and  can  be  synchronized  to  the  same  clock.  Trigger  signals  can  be  distributed  between 
modules.  The  acquisition  of  data  was  simplified  by  assigning  separate  triggers  to  the 
front  and  back  detector  strips.  Data  corresponding  to  coincidence  between  modules  is 
determined  in  the  image  processing  code.  Figure  16  shows  the  clock  and  charge 
collection  strip  distribution  used  for  this  work. 


Figure  16;  Clock  and  strip  distribution 
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The  DGF-4C  modules  have  four  building  blocks:  the  analog  signal  conditioning, 
the  real-time  processing  units  (RTPU),  the  digital  signal  processors  (DSP),  and  the 
CAMAC  interface. 

Analog  Signal  Conditioning:  Each  of  the  analog  inputs  has  its  own  signal 
conditioning  unit,  which  adapts  the  input  signals  to  the  input  voltage  range  of  the  ADC. 
There  is  also  a  software  controllable  gain  setting  at  this  stage.  The  ADC  is  a  waveform 
digitizer. 

RTPU:  There  is  one  RTPU  for  each  channel,  and  it  consists  of  a  field 
programmable-  gate-array  (FPGA)  and  a  first- in  first-out  (FIFO)  memory.  When  data 
arrive  from  the  ADC,  the  RTPU  applies  a  digital  trapezoidal  filter  which  acts  as  a  shaping 
amplifier.  In  addition  to  the  pulse  shaping,  the  RTPU  contains  a  pileup  inspector.  If  the 
logic  detects  a  second  pulse  too  close  to  the  one  being  processed,  so  that  the  pulse  height 
measurement  of  the  first  pulse  would  be  corrupted,  both  of  the  pulses  are  rejected.  If  the 
pulse  is  not  rejected  as  a  pile  up  event,  a  trigger  is  issued  to  notify  the  digital  signal 
processor  (DSP)  that  there  is  raw  data  ready  for  processing.  When  the  trigger  is  issued  to 
the  DSP,  the  RTPU  is  latched  until  it  has  been  serviced  by  the  DSP. 

DSP:  The  DSP  controls  the  operation  of  the  pulse  processing  electronics,  reads 
the  raw  data  from  the  RTPU,  determines  the  pulse  heights,  time  stamps  each  pulse,  and 
prepares  the  data  for  output  to  the  controlling  computer. 

CAMAC  Interface:  The  CAMAC  interface  is  the  means  of  communication  with 
the  DGF-4C. 
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Input  logic  unit  [11] 


The  input  logic  unit  used  is  an  ORTEC  Quad  4-Input  Logic  Unit  Revision  D, 
model  C04020  and  serial  number  382.  The  logic  functions  it  can  perform  are: 
coincidence  (AND),  anti-coincidence  (veto),  fan-in  (OR),  fan-out,  fast  negative  NIM-to- 
TTL  conversion,  and  pulse  lengthening.  For  this  work,  the  fan-in  (OR)  function  was 
utilized.  Each  of  the  four  inputs  accepts  NIM  fast  negative  logic  pulses.  Front-panel, 
three-position  slide  switches  select  the  logic  requirements  separately  for  each  input. 

For  this  work,  the  Busy  Out  ports  on  the  DGF-4Cs  are  each  connected  to  an  input 
on  the  logic  module.  The  logic  module  output  is  distributed  to  each  DGF-4C  Sync  In 
port.  This  setup  allows  both  modules  to  start  and  stop  acquisition  at  the  same  time. 
Coupling  this  with  a  global  clock  for  both  modules  allows  the  determination  of 
coincident  events. 


Power  Distribution  Module  [16] 

The  PDM  used  for  this  effort  is  an  XIA  CAMAC  module  Revision  3  and  serial 
number  108.  The  PDM  provides  +  12  V  and  +  24  V  preamplifier  supply  voltages.  It  is 
capable  of  supplying  power  for  up  to  20  preamplifiers. 


Crate  and  Crate  Controller  [5]  [16] 

The  CAMAC  crate  used  in  this  research  is  a  Wiener  type  UEC  01  VH  /  Typ  342, 
part  number  0305.0102,  and  serial  number  3698023.  It  is  a  25-slot  CAMAC  crate.  The 
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crate  controller  is  a  Jorway  model  13A-2,  serial  number  662.  It  is  the  interfaee  between 
the  CAMAC  crate  and  the  host  eomputer.  The  interfaee  is  via  a  SCSI  eonneetion. 


Collimator 

Proper  eollimation  of  the  incoming  photons  is  crucial  for  producing  an  image.  In 
order  to  obtain  true  spatial  discrimination,  the  fields  of  view  for  the  collimator  holes 
should  not  overlap  at  either  the  souree  or  the  deteetor.  Several  collimators  were  made  in- 
house  for  this  work.  All  were  made  from  AIM-70  (CerroBEND)  using  an  in-house  mold 
design.  AIM  70  is  made  of  50%  bismuth,  26.7%  lead,  13.3%  tin,  and  10%  cadmium.  It 
has  a  melting  point  of  70“  C  whieh  make  it  ideal  for  building  collimators  in-house. 

Several  single  hole  and  three  hole  eollimators  were  used  to  eharaeterize  the 
detector  response.  A  slot  hole  collimator  50  mm  by  50  mm  by  100  mm  with  a  1  mm  slot 
was  used  to  determine  the  deteetor  response  near  the  edges  of  the  charge  colleetion  strips. 
The  characterization  of  the  response  near  the  strip  edges  helped  ensure  the  design  of  the 
final  eollimator  would  provide  proper  spatial  diserimination.  A  parallel  hole  eollimator 
developed  by  Nuclear  Fields  Ine.  for  previous  work  with  the  detector  was  used  to  test  the 
image  processing  eode  before  the  final  eollimator  was  completed.  The  parallel  hole 
collimator  is  made  of  PbSb  alloy  and  consists  of  25  round  holes  3.6  mm  in  diameter.  It  is 
30  mm  by  30  mm  by  50  mm. 

The  eollimator  used  for  final  imaging  was  designed  using  a  Mathematiea  code 
developed  in  house.  The  eode  and  its  output  are  presented  in  Appendix  A.  Collimator 
Design  Code.  The  imaging  collimator  measures  60  mm  by  60  mm  by  190  mm.  It  has  a  5 
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by  5  array  of  square  holes,  2.9  mm  wide  along  their  entire  length.  The  holes  converge 
outward  from  the  detector  to  permit  imaging  objects  larger  than  the  detector  crystal.  The 
hole  configuration  allows  a  5  cm  radius  source  to  be  imaged  from  a  distance  of  7 1  cm 
from  the  source  side  of  the  collimator.  This  configuration  was  chosen  to  simulate 
imaging  a  nuclear  weapon  pit  while  it  is  in  a  storage  cask.  Figure  17  shows  the 
configuration  used  to  properly  align  the  converging  hole  collimator  with  the  detector  for 
imaging  operations. 


Figure  17:  Imaging  configuration  using  converging  hole  collimator 


Image  Processing  Code 

In  order  to  produce  an  image,  the  raw  data  collected  from  the  DGFs  was 
processed  using  a  processing  code  written  in  the  MATLAB  computer  language.  The 
code  contains  9  main  algorithms  for  processing  the  raw  data.  The  9  steps  are  listed 
below,  and  the  full  code  is  presented  in  Appendix  B.  Image  Processing  Code. 
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1.  Read  in  the  raw  data  and  user  inputs.  The  raw  data  is  saved  using  the  IGOR 
software  as  a  text  fde.  The  user  is  prompted  to  input  up  to  two  energy  windows  and 
ealibration  for  the  image.  A  default  ealibration  is  also  available.  The  user  also  has  the 
option  to  use  three  strip  events  in  addition  to  two  strip  events. 

2.  Sort  the  data.  An  event  data  set  eonsists  of  the  hit  pattern,  timestamp  and 
ehannel  number  eorresponding  to  the  energy.  Eaeh  event  is  sorted  into  one  of  four 
matrices  depending  on  which  module  it  is  for  and  whether  it  is  a  single  strip  or  two  strip 
coincident  event. 

3.  Each  single  rear  strip  event  is  then  compared  to  each  single  front  strip  event  to 
determine  events  coincident  in  both  time  and  energy. 

4.  When  coincident  events  are  found  in  step  3,  each  event  hit  pattern  is  used  to 
determine  which  strip  received  the  event.  An  event  hit  is  then  assigned  to  the  appropriate 
image  pixel  by  adding  one  to  its  value. 

5.  Two  strip  events  for  the  rear  strips  are  compared  to  each  single  front  strip  event 
to  determine  three  strip  events  coincident  in  both  time  and  energy. 

6.  When  three  strip  coincident  events  are  found  in  step  5,  each  event  hit  pattern  is 
used  to  determine  which  strips  received  the  event.  An  event  hit  is  then  assigned  to  the 
image  pixel  which  recorded  its  portion  of  the  charge  first  by  adding  one  to  its  value. 

7.  Two  strip  events  for  the  front  strips  are  compared  to  each  single  rear  strip  event 
to  determine  three  strip  events  coincident  in  both  time  and  energy. 

8.  When  three  strip  coincident  events  are  found  in  step  7,  each  event  hit  pattern  is 
used  to  determine  which  strips  received  the  event.  An  event  hit  is  then  assigned  to  the 
image  pixel  which  recorded  its  portion  of  the  charge  first  by  adding  one  to  its  value. 
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9.  The  final  image,  as  seen  by  the  deteetor,  is  then  shown  along  with  a  surfaee 
plot  of  the  image  and  a  histogram  showing  the  number  of  events  reeorded  for  eaeh  energy 
throughout  the  speeified  energy  windows. 

In  steps  6  and  8  above,  the  assignment  of  a  three  strip  event  to  a  single  pixel  is 
eontrolled  by  the  timing  of  the  pulses  produeed  on  adjaeent  strips.  The  timing  resolution 
of  the  DGF  modules  is  25  ns.  This  means  that  pulses  on  adjaeent  strips  whieh  reeord  the 
same  time  must  have  oeeurred  within  25  ns  of  eaeh  other.  Therefore,  adjaeent  strip 
events  whieh  oeeur  at  the  same  reeorded  time  eannot  be  used  with  this  system.  However, 
if  the  time  between  adjaeent  strip  pulses  is  more  than  25  ns  then  it  is  assumed  that  the 
gamma  ray  entered  the  deteetor  at  the  strip  whieh  reeorded  its  pulse  first.  This 
assumption  is  possible  beeause  the  timing  trigger  for  the  pulses  is  based  off  a  threshold 
value  and  not  a  pereentage  of  the  individual  pulse. 


Radiation  Sources 

Activated  Tin  Source  [14].  In  order  to  simulate  Pu-239,  a  distributed  radiation 
souree  was  needed  that  had  photons  in  the  375  and  413  keV  ranges.  The  easiest  souree  to 
get  that  satisfied  the  neeessary  requirement  was  neutron  aetivated  tin.  The  5  Sn  strips 
were  5  mm  wide  and  25  mm  long  and  approximately  1 .7  mm  thiek.  The  five  Sn  strips 
were  aetivated  at  the  Ohio  State  University  researeh  reaetor  on  28  August  2003.  The 
strips  had  masses  ranging  from  1 .493g  to  1 .5 12g.  All  five  strips  were  plaeed  in  the  OSU 
reaetor  eentral  irradiation  faeility  (GIF)  for  two  and  a  half  hours,  while  the  reaetor  was 
run  at  90%  power.  After  the  irradiation  of  the  sourees  was  eompleted,  the  reaetor  was 
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shut  down  and  the  sources  remained  in  the  CIF  overnight  to  allow  the  extremely  short 
lived  isotopes  to  decay.  Then  the  sources  were  removed  from  the  CIF  and  were  placed 
into  a  pig,  where  they  were  allowed  to  sit  for  three  weeks  to  let  the  short-lived  activation 
products  decay  away.  After  this  cooling  off  period,  only  long-lived  activation  products 
remained.  The  long-lived  activation  products  resulting  from  neutron  activated  Sn  are 
given  in  Table  1.  The  activities  of  the  strips  are  unknown.  Because  of  its  relatively  short 
half-life,  Sn-1 17m  was  not  detected  during  this  work.  However,  the  relatively  short  lived 
Sb-124  was  still  present  because  there  was  such  a  large  quantity  present  following  the 
neutron  activation.  A  spectrum  of  one  neutron-activated  Sn  sample,  as  seen  in  the  DGF 
controlling  software,  recorded  by  a  single  strip  from  the  HPGeDSSD  is  shown  in  Figure 
18.  The  collection  time  for  the  spectrum  was  2  hours.  The  spectrum  was  collected  on  3 
December  2004  (463  days  after  activation). 

Table  1;  Activated  Sn  strip  gamma  emissions 


Nuclide 

Gamma-ray 
Energy  [keV] 

Branching 
Ratios  (%) 

Half-life 

Sn  117m 

158.91 

86.0 

14  days 

Sn  113 

255.30 

1.8 

115  days 

391.73 

64.97 

Sb  124 

602.7 

98.26 

60  days 

Sb  125 

176.57 

6.82 

2.76  years 

427.89 

30.0 

463.34 

10.5 

600.56 

17.9 

606.34 

5.03 

635.88 

11.31 
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Figure  18;  Activated  Sn  strip  spectrum 

Simulated  Weapon  Pit  Source  [9]:  In  order  to  test  the  ability  of  the  system  to 
image  a  weapon  pit  in  a  configuration  similar  to  that  encountered  during  on  site  weapon 
inspections,  a  distributed  Cs-137  bowl  source  was  used.  The  source  was  designed  and 
specified  by  Mr.  Mike  Nelson,  a  previous  AFIT  student.  The  bowl  was  constructed  from 
a  polymer  material  to  form  a  hemisphere  shell.  The  outer  and  inner  surfaces  are  0.3 
inches  (7.62  mm)  thick  surrounding  a  hollow  cavity  constructed  to  be  0.3  inches  (7.62 
mm)  across.  The  outer  radius  of  the  hollow  cavity  is  5  cm.  The  configuration  produces  a 
5  cm  radius  hemispherical  source  cavity  which  is  7.62  mm  thick.  The  interior  cavity 
contains  a  Cs-137  source  suspended  in  approximately  54  ml  of  IM  HCL.  The  reference 
activity  of  the  source  was  650  //Ci  on  10  January  2003. 

Shortly  after  fabrication,  the  source  was  found  to  be  leaking  the  Cs-137  solution. 
This  led  to  the  source  being  sealed  in  a  plastic  bag  for  safe  handling  and  storage.  I  was 


37 


unable  to  determine  the  amount  of  souree  solution  whieh  had  leaked  out,  but  it  appeared 
from  the  aequired  images  that  enough  was  lost  to  affeet  the  symmetry  of  the  souree. 

Standard  Sources:  The  following  standard  sources  were  used  for  determining 
the  calibration,  efficiency  and  resolution  of  the  detector.  The  sources  were  also  used  for 
aligning  the  collimator  with  the  detector  and  showing  the  ability  to  image  various  sources 
at  the  same  time  using  multiple  energy  windows. 

Table  2:  Radiation  source  used 


Source 

(Reference  date) 

Nuclide 

Gamma-ray 

Energy 

fkeV] 

Initial 

Activity 

(hCi) 

Branching 
Ratio  (%) 

Half-life 

Multi  Nuclide  - 
T108  Isotope 
Products 
Laboratories 
(15  Jul  04) 

Am-241 

60 

0.02912 

36.0 

432.17 

Cd-109 

88 

0.4260 

3.63 

462.6  days 

Co-57 

122 

0.01529 

85.6 

271.79  days 

Ce-139 

166 

0.01963 

79.9 

137.64  days 

Hg-203 

279 

0.06120 

81.5 

46.595  days 

Sn-113 

392 

0.07482 

64.9 

115.09  days 

Sr-85 

514 

0.09425 

98.4 

64.849  days 

Cs-137 

662 

0.06875 

85.1 

30.17  years 

Y-88 

898 

0.1481 

94.0 

106.63  days 

Co-60 

1173 

0.08091 

99.86 

5.272  years 

Co-60 

1333 

0.08091 

99.98 

5.272  years 

Y-88 

1836 

0.1481 

99.4 

106.63  days 

Na-22-T107 
Isotope  Products 
Laboratories 
(15  Feb  03) 

Na-22 

511 

0.8796 

99.94 

950.8  days 

1275 

0.8796 

99.94 

950.8  days 

Eu-152-TllO 
Isotope  Products 
Laboratories 
(15  Dec  04) 

Eu-152 

244 

10.14 

7.5 

4933  days 

344 

10.14 

26.6 

4933  days 

CS-I37-T089 
Isotope  Products 
Laboratories 
(I  Aug  98) 

Cs-137 

662 

10.02 

85.1 

30.17  years 
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Controlling  Software  /  SCSI  Interface 


The  controlling  software  used  for  the  DGF-4C  modules  was  the  XIA  DGF-4C 
viewer,  release  number  3.04.  This  software  runs  inside  the  IGOR  Pro  software  which 
was  supplied  by  XIA  with  the  purchase  of  the  DGF-4C  modules.  The  controlling 
software  provides  the  user  a  platform  from  which  all  settings  can  be  adjusted  for  the 
entire  pulse  processing  system.  It  also  provides  an  environment  for  analyzing  the 
resulting  data.  The  SCSI  interface  card  used  in  the  PC  is  an  Adaptec  2930  SCSI  card. 
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IV.  Procedure 


Chapter  Overview 

In  order  to  help  optimize  the  imaging  system,  a  eollimator  design  eode  was 
developed.  This  allowed  for  design  of  several  eollimators  in  order  to  test  the  deteetor 
response  under  variable  eonfigurations.  An  image  proeessing  eode  was  also  developed 
early  on.  Both  eodes  were  updated  throughout  the  experiment  to  meet  the  required  needs 
as  they  were  identified. 

Prior  to  any  eollimator  testing,  the  equipment  was  set  up  and  the  proper  hardware 
and  software  settings  required  to  ereate  an  image  were  determined.  This  ineluded 
charaeterizing  the  deteetor  by  determining  the  noise  present,  the  proper  energy  ealibration 
and  effieiency  of  the  deteetor.  With  eaeh  strip  calibrated  to  the  same  energy  scale, 
several  test  collimators  were  used  to  determine  the  optimum  size  parameters  of  the  final 
collimator.  This  work  also  included  investigating  the  possibility  of  using  holes  located  at 
the  edge  of  each  pixel.  Once  a  feasible  collimator  design  was  developed,  mold  drawings 
were  made  and  the  mold  was  fabricated.  The  collimator  was  then  built  and  put  into  the 
system  for  imaging. 


Collimator  Design  and  Fabrication 

The  Mathematica  code  developed  for  the  collimator  design  is  shown  in  Appendix 
A.  The  code  allows  the  user  to  optimize  the  design  by  varying  the  size  of  the  holes, 
collimator  length,  and  focal  length  of  the  converging  holes.  All  holes  are  required  to  be 
the  same  size  for  the  entire  length  of  the  collimator. 
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Theoretically,  a  converging  collimator  using  edge  pixels  could  be  used  to  increase 
the  spatial  resolution  of  the  source  size  and  shape  while  increasing  the  size  of  the  source 


as  compared  to  the  previous  work  with  this  detector.  One  such  example  is  the  partial 
nine  by  nine  array  of  65  holes  shown  in  Figure  19.  Black  filled  circles  represent  the 
source  face  configuration  and  the  white  open  circles  represent  the  detector  face 
configuration.  Figure  20  shows  the  resulting  fields  of  view  for  each  collimator  hole  in 
the  center  of  the  detector  crystal. 


Figure  19;  Proposed  collimator  hole  pattern 


Figure  20:  Proposed  hole  pattern  fields  of  view  on  detector  center  plane 
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A  test  collimator,  shown  in  Figure  21,  was  built  to  test  the  feasibility  of  using 
edge  pixel  collimator  holes  for  gamma  imaging.  This  was  done  by  using  a  tin  strip 
source  collimated  with  a  1  mm  wide  gap  along  the  entire  25  mm  length  of  the  strip.  The 
response  of  strips  F2  and  F3  was  investigated  with  the  collimated  source  located  at 
various  locations  near  the  edge  between  the  strips.  I  expected  to  be  able  to  find  a  range 
near  the  gap  in  which  each  strip  recorded  very  few  full  energy  peaks.  This  would  then 
determine  the  size  limit  of  an  edge  pixel  hole.  If  the  edge  pixel  records  a  full  energy  peak 
on  a  single  strip  then  it  would  be  assigned  to  the  center  strip  pixel  for  that  strip  and 
therefore  not  provide  proper  spatial  resolution  of  the  source.  As  discussed  in  Chapter  5, 
the  recorded  data  showed  that  the  edge  pixel  holes  would  have  to  be  too  small  to  provide 
enough  efficiency  for  imaging.  This  led  me  to  use  only  center  pixel  holes  for  the  final 
collimator  design. 


Figure  21 :  Edge  pixel  test  collimator 
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The  original  collimator  design  code  utilized  round  holes.  However,  analysis  of 
the  design  showed  that  using  square  holes  would  increase  the  field  of  view  areas  on  the 
source  and  the  detector  by  approximately  25%.  This  increase  in  the  volume  of  the 
detector  utilized  and  the  volume  of  the  source  visible  to  the  detector  results  in  shorter 
data  acquisition  times  and  a  higher  efficiency  of  the  system. 

The  primary  objective  for  the  collimator  design  was  to  maximize  the  amount  of 
the  detector  used  while  also  maximizing  the  amount  of  the  source  seen.  In  addition, 
proper  spatial  resolution  must  be  achieved  by  ensuring  the  fields  of  view  within  the 
detector  and  at  the  source  do  not  overlap.  Figure  22  shows  the  optimized  collimator 
design  hole  pattern.  The  collimator  is  190  mm  long  with  2.9  mm  square  holes  and  a  1.16 
m  focal  length.  The  black  holes  represent  the  hole  configuration  on  the  side  nearest  to 
the  detector.  The  gray  holes  represent  the  hole  configuration  on  the  side  nearest  to  the 
source.  Appendix  A  presents  the  design  code  and  its  output  for  the  optimized  collimator 
configuration. 


Figure  22:  Optimized  collimator  hole  pattern 

Figure  23  shows  the  resulting  fields  of  view  superimposed  on  the  detector  pixel 
array  and  simulated  weapon  pit  source.  In  the  left  image,  the  black  areas  represent  the 
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fields  of  view  on  the  rear  face  of  the  detector.  The  gray  areas  represent  the  fields  of  view 
on  the  front  face  of  the  detector.  In  the  right  image,  the  circular  area  represents  the  cross 
section  of  a  5  cm  radius  source  centered  710  mm  from  the  front  face  of  the  collimator. 

As  Figure  23  shows,  each  collimator  hole  has  a  field  of  view  measuring  approximately 
2.5  mm  on  each  side.  Thus  the  origin  of  each  photon  emitted  can  be  known  to  within 
approximately  6.25  cm  . 


Figure  23:  Collimator  fields  of  view  on  the  detector  and  source 


Once  the  desired  collimator  dimensions  were  determined,  initial  mold  design 
drawings  were  developed.  These  drawings  were  then  taken  to  a  local  machine  shop, 
Superior  Machining  Inc.,  for  final  mold  design  and  initial  fabrication.  Final  mold 
fabrication  included  development  of  the  25  square  rods  used  to  produce  the  holes  and 
cutting  out  the  square  holes  in  the  endplates  for  proper  placement  of  the  rods.  These  final 
fabrication  steps  were  complete  by  the  AFIT  Model  Fabrication  Shop.  Figure  24  shows 
the  components  of  the  collimator  mold  assembly. 
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Eigure  24;  Collimator  mold  assembly 


Upon  completion  of  the  mold,  the  final  collimator  was  poured  in  house  using 
CerroBEND.  The  first  step  for  pouring  the  collimator  was  to  assemble  the  mold.  The 
detector  and  source  endplates  were  attached  to  the  trough  and  the  25  rods  were  placed  in 
the  endplates.  Due  to  limitations  in  the  electrical  discharge  machining  (EDM)  device 
required  to  cut  the  holes  in  the  endplates,  the  square  holes  had  to  be  cut  perpendicular  to 
the  endplate  faces.  This  led  to  the  requirement  of  oversized  holes  which  allow  each  rod 
to  be  placed  at  the  proper  angle.  The  use  of  oversized  holes  produced  gaps  on  some 
edges  between  the  endplates  and  the  rods.  In  order  to  keep  the  CerroBEND  from  leaking 
out  the  gaps,  I  packed  the  gaps  with  white  lithium  grease.  Each  rod  was  also  coated  with 
a  thin  film  of  white  lithium  grease  to  aid  in  removal  after  the  CerroBEND  solidified. 
Eigure  25  shows  the  collimator  mold  fully  assembled. 
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Figure  25:  Fully  assembled  eollimator  mold 

The  CerroBEND  used  to  make  the  eollimator  melts  at  70°  C  so  I  melted  it  at  85° 

C  to  ensure  I  would  have  time  to  pour  it  before  it  began  to  solidify.  Onee  the  material 
was  melted,  I  then  heated  the  assembled  mold  in  the  oven  for  1 5  minutes  in  order  to  keep 
the  material  from  solidifying  too  fast  as  it  hit  the  mold.  The  mold  was  then  removed 
from  the  oven  and  the  melted  CerroBEND  was  poured  in.  The  full  mold  was  plaeed  back 
into  the  oven  to  help  remove  some  of  the  air  bubbles  in  the  CerroBEND  formed  while 
pouring.  The  mold  was  then  removed  from  the  oven  and  the  cooling  plate  was  placed  on 
top  of  the  mold  to  aid  in  uniform  cooling.  A  small  fish  tank  air  pump  was  placed  in  the 
pan  next  to  the  mold  as  it  cooled  to  provide  vibrations  which  aid  in  removing  air  bubbles. 
As  the  assembly  cooled,  I  tapped  the  25  rods  back  and  forth  to  help  keep  the  CerroBEND 
from  adhering  to  the  rods. 

Once  the  CerroBEND  had  cooled  sufficiently,  I  then  had  to  remove  the  25  rods. 
Unlike  most  metals,  as  CerroBEND  cools  it  actually  expands  rather  than  contracts.  This 
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led  to  some  difficulty  in  removing  the  rods,  even  with  the  application  of  the  white  lithium 
grease.  The  rods  had  to  be  driven  out  using  a  small  hammer  and  rods  of  square  key  stock 
steel  slightly  smaller  than  the  dimensions  of  the  rods.  Once  the  rods  were  removed  and 
the  holes  were  cleaned  out,  the  collimator  was  placed  into  a  base  assembly  fabricated  to 
allow  it  to  be  attached  to  a  translation  stage  for  proper  alignment  with  the  detector. 

Equipment  Setup 

The  first  step  in  setting  up  the  equipment  was  to  ensure  that  all  of  the  internal 
switch  settings  and  jumper  settings  for  the  Jorway  crate  controller  were  in  the  correct 
position.  There  are  four  internal  switches  that  were  all  placed  in  the  OFF  position,  a 
jumper  was  placed  spanning  X2  and  X3,  and  it  was  verified  that  the  controller  was 
internally  terminated.  Then  the  controller  was  placed  into  the  two  right-most  slots  of  the 
crate,  and  the  PDM  was  placed  into  the  two  left-most  slots. 

Initially,  only  a  single  DGF-4C  module  was  used  to  verify  that  the 
communication  with  the  host  computer  was  functioning  properly.  All  internal  jumper 
settings  were  verified  as  being  positioned  correctly  for  a  single  module.  Once  the  ver 
3.04  XIA  software  was  installed  on  the  host  computer,  it  was  connected  to  the  controller 
with  a  SCSI  cable  and  the  module  was  successfully  booted.  During  the  start  up  process 
the  SCSI  bus  ID  number  and  crate  number  are  required.  The  bus  ID  number  is  2  and  the 
crate  number  is  set  to  3  the  way  the  system  is  currently  configured. 

Once  the  module  was  booted  up  a  signal  pulse  generator  was  used  as  a  signal 
input  to  become  familiar  with  the  controlling  software.  Figure  26  shows  the  setup  used 
during  this  familiarization. 
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Figure  26:  Experiment  test  setup 


After  becoming  familiar  with  the  software,  two  DGF-4C  modules  were  placed 
into  the  crate  and  connected  across  their  backplane  to  distribute  the  master  clock,  and  the 
internal  jumpers  were  all  placed  in  the  XIA  recommended  positions.  The  following 
notation  is  used  throughout  this  project,  module  1  (rear  strips)  is  the  right  XIA  module 
and  module  2  (front  strips)  is  the  left  XIA  module.  The  following  jumper  settings  were 
used  when  both  modules  were  placed  in  the  crate  together. 

Table  3:  DGF-4C  module  jumper  settings 


Module  1  -  1244 

Module  2  -  1240 

Jumper  1 

Installed 

Removed 

Jumper  2 

Installed 

Removed 

Jumper  3 

Installed 

Installed 

Jumper  4 

Installed 

Installed 

Jumper  5 

Board  Clock 

External 
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After  the  modules  were  placed  in  the  crate,  a  ribbon  cable  provided  by  XIA  was 
used  to  connect  the  modules  together  via  the  8  pin  connectors  on  the  back  of  each 
module.  Figure  27  shows  the  backplane  connections  as  viewed  from  behind  the  modules. 


1244  1240 


Figure  27:  DGF-4C  module  backplane  connection 

After  the  modules  had  been  hooked  together  across  the  backplane,  the  logic 
module  was  connected  to  the  front  of  the  modules.  The  DGF  ‘busy  out’  connectors  were 
hooked  to  the  input  of  the  logic  module,  and  the  DGF  ‘synch  in’  connectors  were  hooked 
to  the  output  of  the  logic  module.  Once  all  connections  had  been  made,  the  two  modules 
were  booted  up  and  four  pulse  generators  were  connected  to  the  inputs  of  the  DGF-4C’s. 
All  four  of  the  signals  were  recorded  by  the  DGF  software  in  MCA  mode,  showing  that 
the  system  was  functioning  properly  and  ready  for  the  detector  outputs  to  be  connected  to 
the  DGF-4C  modules. 

The  dewar  for  the  HPGeDSSD  was  verified  to  be  full  and  the  detector  was 
connected  to  an  ORTEC  659  high  voltage  power  supply.  The  power  supply  internal 
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board  was  placed  in  the  positive  polarity  position,  and  in  a  position  to  allow  automatie 
shut  off.  The  high  voltage  shut  down  (HVSD)  eonneetor,  on  the  deteetor,  was  hooked  to 
the  baek  of  the  power  supply  at  its  automatie  shut  down  eonneetion.  This  is  a  safety 
feature  whieh  will  foree  a  shutdown  of  the  high  voltage  power  supply  if  the  deteetor 
temperature  is  allowed  to  rise  above  liquid  nitrogen  levels. 

The  high  voltage  supply  line  was  hooked  to  the  HV-In  eonneetor  on  the  deteetor, 
and  the  PDM  was  eonneeted  to  the  preamplifier  supply  eonneetion  on  the  deteetor.  The 
deteetor  outputs  were  eonneeted  to  the  DGF-4C  modules.  For  this  work,  strips  F1-F4  and 
R1-R4  were  used  to  produee  the  desired  four  by  four  array.  Then  the  high  voltage  power 
was  turned  on  and  slowly  inereased  to  +1000  volts  in  100  volt  inerements.  At  eaeh 
inerement  of  100  volts,  the  deteetor  outputs  were  observed  using  the  oseilloseope  feature 
in  the  DGF  software  to  ensure  that  as  power  was  inereased  the  noise  on  eaeh  ehannel 
deereased.  Table  4  lists  whieh  eonneetions  were  made  between  deteetor  outputs  and 
DGF-4C  module  inputs. 

Table  4:  Deteetor  output  to  DGF-4C  module  eonneetions 


Deteetor 

Module 

Channel 

Output 

Number 

Number 

FI 

2 

0 

F2 

2 

1 

F3 

2 

2 

F4 

2 

3 

R1 

1 

0 

R2 

1 

1 

R3 

1 

2 

R4 

1 

3 
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The  XIA  notation  for  connectors  is:  the  top  connector  on  the  DGF-4C  modules  is  channel 
0,  and  the  bottom  connector  is  channel  3. 


Noise 

Noise  present  on  each  detector  strip  can  be  viewed  using  the  oscilloscope  function 
within  the  DGF  software.  During  application  of  the  high  voltage  bias  to  the  detector,  the 
noise  can  be  observed  to  ensure  it  decreases  as  the  bias  increases.  Failure  of  the  noise  to 
decrease  would  indicate  a  problem  with  the  bias  being  applied. 


Charge  Sharing  Between  Strips 

To  ensure  complete  charge  collection  was  being  achieved,  all  ten  detector  outputs 
were  viewed  simultaneously  on  oscilloscopes.  Two  four  channel  and  one  two  channel 
oscilloscopes  were  used.  All  ten  of  the  channels  were  set  to  trigger  off  of  the  channel 
displaying  the  detector  output  FI .  Each  pulse  could  then  be  matched  to  a  single  pulse  or 
multiple  pulses  occurring  on  the  opposite  face  of  the  detector. 


Software  Settings 

Preamplifier  decay  constants:  One  of  the  user  inputs  within  the  XIA  software  is 
the  preamplifier  decay  constant.  This  is  used  by  the  software  to  correct  for  pulses  that 
arrive  before  the  previous  pulse  has  returned  to  baseline.  Since  the  response  of  each  strip 
is  slightly  different,  a  preamplifier  constant  is  required  for  each  strip.  Optimization  of  the 
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preamplifier  deeay  eonstant  is  performed  using  the  DGF  software.  Ten  values  of  the 
eonstant  are  ealeulated  and  the  average  value  is  shown  to  the  user  for  aeoeptance  or 
rejection. 

Gain  Settings:  Proper  energy  calibration  of  each  strip  is  crucial  when  imaging  in 
a  narrow  energy  range.  This  requires  that  each  strip  be  carefully  calibrated  to  the  same 
energy  scale.  This  is  accomplished  by  adjusting  each  strip  gain  setting  so  that  it  will 
place  the  centroid  of  a  full  energy  peak  (FEP)  at  a  specified  channel.  The  FEP  chosen  for 
the  energy  calibration  was  the  392  keV  peak  from  Sn-1 13.  The  location  for  this  peak  was 
chosen  to  be  channel  20200. 

Energy  Filter  Rise  Time:  Optimization  of  energy  resolution  was  performed  by 
determining  the  energy  filter  rise  times  which  resulted  in  a  minimized  EWHM  for  each 
strip.  This  was  done  at  the  suggestion  of  the  XIA  users  guide  for  maximizing 
performance.  Eourteen  data  points  were  taken  at  different  energy  filter  rise  times,  and 
each  spectra  resolution  was  measured.  All  measurements  were  taken  using  a  Cs-137 
source  with  a  collection  time  of  600  seconds. 

Common  XIA  Software  Settings:  Each  of  the  DGE-4C  modules  and  channels 
has  other  settings  that  have  to  be  set  within  the  software.  Many  of  the  settings  are  the 
same  for  all  of  the  channels.  The  settings  that  are  constant  for  each  channel  are  listed  in 
Table  5. 
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Table  5:  Common  XIA  software  settings 


Trigger  filter  rise  time  (ps) 

0.1 

Trigger  fdter  flat  top  (ps) 

0.1 

Energy  filter  flat  top  (ps) 

2.4 

Threshold 

200 

Cutoff 

5000 

Binning  factor 

1 

Coincidence  pattern 

115E 

Resolution  Determination 

As  deseribed  above,  the  energy  filter  rise  time  was  optimized  for  eaeh  strip  in  an 
attempt  to  optimize  the  resolution  of  the  deteetor  as  a  whole.  For  the  purpose  of  imaging, 
the  resolution  was  measured  at  391.7  keV  with  emissions  from  the  Sn-1 13  isotope,  in  the 
activated  Sn  source,  with  a  counting  time  of  36,000  seconds.  The  resolution  of  each  of 
the  charge  collection  strips  was  measured  using  a  normal  distribution  fit  that  is  performed 
by  the  XIA  software.  The  fit  was  used  because  it  automatically  subtracts  out  the 
background.  The  fit  gives  both  the  absolute  FWHM,  and  the  percentage  FWHM.  For 
purposes  of  reporting  the  resolution,  the  absolute  FWHM  value  is  used. 


Efficiency  Calculation 

To  determine  the  efficiency  of  the  detector,  four  radiation  sources  of  known 
activities  were  used,  for  a  total  of  six  full  energy  peaks.  The  efficiency  decreases  with 
increasing  photon  energy.  Depending  on  its  activity,  each  source,  one  at  a  time,  was 
placed  at  a  specified  distance  along  the  axis  of  the  detector.  The  counting  time  for  each 
source  was  also  varied  depending  upon  the  source  activity.  For  the  calculation  of  the 
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solid  angle,  I  assumed  a  right  eircular  detector  with  a  radius  that  had  a  surface  area  equal 
to  the  surface  area  of  the  ORTEC  detector.  I  then  established  a  setup  with  a  source-to- 
detector  separation  much  larger  than  the  source  dimension.  This  allowed  me  to  use  the 


following  equation  for  calculating  the  solid  angle  [6] 


f 

Q.  =  ln  1 

V 


d 


4d 


J 


(3) 


where: 

n=  Solid  angle 

d  =  Distance  between  the  source  and  detector  [cm] 
a  =  Detector  radius  [cm] 


The  intrinsic  efficiency  was  then  calculated  using  the  following  equation  [6] 


counts  recorded 


A-t-b- 


An 


(4) 


where: 

Sint=  Intrinsic  efficiency 
A  =  Activity  of  the  source  (Bq) 
t  =  Counting  time  (sec) 
h=Branching  ratio 
Q.  =  Solid  angle 


Collimator  Alignment 

Proper  alignment  of  the  converging  collimator  with  the  detector  strips  is  required 
to  create  an  image  with  proper  spatial  resolution.  Proper  alignment  is  defined  as  aligning 
each  collimator  hole  with  only  one  corresponding  detector  pixel.  This  ensures  that  a 
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gamma  photon  producing  an  event  in  a  particular  pixel  originated  from  within  a  known 
volume  of  the  source.  Initially,  I  assumed  that  the  detector  was  located  in  the  eenter  of 
the  eylindrieal  eryostat  housing.  The  detector  manufacturer  reports  that  the  front  deteetor 
strips  are  between  10  and  15  mm  from  the  front  of  the  cryostat  housing.  This  uncertainty 
required  that  proper  eollimator  alignment  be  performed  experimentally. 

The  collimator  was  placed  into  a  base  assembly  and  seeured  to  a  three  axis 
adjustable  stage.  The  collimator  was  then  aligned  with  the  approximated  center  of  the 
eryostat  housing,  and  spectra  of  a  point  Eu-152  source  loeated  at  the  eenter  eollimator 
hole  were  taken.  Figure  28  shows  the  eonfiguration  used  for  this  step.  The  actual  Eu-152 
source  is  only  5mm  diameter  so  it  is  only  visible  through  one  collimator  hole  when 
properly  positioned.  The  eollimator  was  then  repositioned  according  to  the  resulting 
spectra.  This  process  was  repeated  until  the  source  peak  was  seen  only  on  the  two  strips 
corresponding  to  the  center  pixel,  namely  strips  R3  and  F3. 


Figure  28:  Center  pixel  alignment  configuration 
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The  next  step  was  to  ensure  proper  alignment  of  the  remaining  pixels.  This  was 
performed  by  plaeing  the  Eu-152  souree  at  one  of  the  eollimator  eorner  holes  and 
aequiring  new  speetra.  The  distanee  from  the  eollimator  to  the  eryostat  housing  was  then 
adjusted  aeeording  to  the  resulting  image.  This  proeess  was  repeated  until  the  image  of 
the  point  souree  was  loeated  in  the  appropriate  pixel.  Proper  alignment  was  then  verified 
by  plaeing  sourees  at  the  opposite  eorner  hole  of  the  eollimator  and  ensuring  that  the 
eollimator  hole  eovered  only  one  deteetor  pixel.  Onee  the  eollimator  alignment  was 
eomplete,  data  aequisition  was  ready  to  begin. 

Imaging 

To  ereate  images,  the  raw  data  were  eolleeted  by  the  deteetor,  and  proeessed  by 
the  MATLAB  image  proeessing  eode  presented  in  Appendix  B.  The  resulting  images 
were  then  analyzed  in  order  to  determine  eertain  eharaeteristies  of  the  system. 

Data  Collection;  For  the  eolleetion  of  imaging  data,  the  XIA  software  must  be 
set  to  take  list  mode  data.  In  the  run  tab  of  the  XIA  software,  the  run  type  must  be  set  to 
energy  and  time  only.  In  this  list  mode,  only  the  hit  patterns,  timestamps  and  energies  are 
stored  for  eaeh  event.  During  operation,  the  memory  in  the  level  1  buffer  for  eaeh  input 
ean  only  store  up  to  744  events,  and  then  it  must  dump  the  data  to  a  larger  memory.  Eaeh 
of  the  eolleetions  of  744  events  is  referred  to  as  a  spill.  Within  the  XIA  software  one  ean 
speeify  a  maximum  time  for  eaeh  spill,  assuming  that  the  buffer  is  not  filled  before  the 
time  runs  out.  One  also  must  speeify  the  number  of  spills  to  be  stored  during  data 
eolleetion.  Onee  the  data  eolleetion  is  eomplete,  the  DGF-4C  module  downloads  all  of 
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the  data  to  the  PC  as  a  binary  file.  Within  the  XIA  software,  the  data  ean  be  eonverted 
from  its  binary  form  and  saved  as  a  text  delimited  file,  ft  is  this  file  that  is  read  in  by  the 
image  proeessing  eode. 

Image  Processing  Code;  Chapter  3  presents  the  key  algorithms  found  in  the 
eode.  The  aetual  eode  is  presented  in  Appendix  B.  Applieation  of  the  eode  requires  a 
few  user  inputs.  First,  the  user  must  ehoose  how  many  energy  windows  to  image  within. 
The  user  must  then  ehoose  to  use  only  two  strip  events  or  both  two  and  three  strip  events. 
Next,  the  user  seleets  the  upper  and  lower  limits  of  the  energy  windows.  The  last  user 
input  is  whether  to  aeeept  the  default  energy  ealibration  or  input  a  new  one.  The  eode 
then  proeeeds  to  read  in  and  proeess  the  raw  data.  The  final  output  of  the  eode  ineludes  a 
figure  window  with  the  four  by  four  pixel  image  of  the  souree  as  seen  from  the  deteetor 
along  with  a  surfaee  plot  of  the  image  and  a  histogram  showing  the  number  of  eounts  for 
eaeh  energy  within  the  speeilied  range.  The  resulting  images  ean  be  exported  to  a  desired 
loeation  in  various  formats  for  further  analysis. 

Image  Analysis:  The  image  proeessing  eode  proved  to  be  a  valuable  tool  for 
analyzing  eertain  features  of  the  images.  The  primary  areas  of  interest  for  the  image 
analysis  ineluded  the  utilization  of  two  strip  events  versus  three  strip  events,  the  ability  to 
distinguish  a  point  souree  from  a  distributed  souree,  and  the  ability  to  image  using 
multiple  energy  windows.  Data  were  taken  for  several  image  eonfigurations  ineluding 
single  souree  and  multiple  souree  setups.  The  images  were  then  proeessed  using  a 
variety  of  energy  windows  in  order  to  explore  the  above  mentioned  areas  of  interest. 
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V.  Results 


Chapter  Overview 

This  chapter  presents  the  results  of  the  deteetor  eharaeterization,  the  statisties 
involved  with  the  imaging,  and  the  aetual  imaging  done  with  the  deteetor. 

Noise 

When  measuring  the  noise  present  on  eaeh  of  the  ehannels,  it  was  found  that  the 
R2  pre-amplifier  output  had  an  approximately  10%  larger  noise  eomponent  than  any  of 
the  other  eharge  eolleetion  strip  outputs.  The  oseilloseope  built  into  the  DGF-4C 
eontrolling  software  was  used  for  these  measurements.  All  ehannels  showed  noise 
eharaeteristie  of  low  frequeney  meehanieal  or  seismie  vibrations.  One  possible 
explanation  to  the  inerease  in  noise  on  ehannel  R2  is  that  the  pre-amplifier  is  more 
sensitive  to  the  above  mentioned  vibrations.  The  extra  noise  on  this  strip  does  not 
degrade  the  energy  resolution  eonsiderably.  Figure  29  shows  a  eomparison  between  the 
eharge  eolleetion  strip  R2  noise  level  and  the  R4  strip  noise  level.  The  other  six  eharge 
eolleetion  strips  have  noises  that  are  eomparable  to  the  R4  strip. 
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Figure  29:  Pre-amplifier  noise  eomparison  between  two  strips 


Charge  Collection/Sharing 

To  ensure  that  full  eharge  eolleetion  is  aeeomplished,  the  preamplifier  outputs 
were  viewed  on  multiple  oseilloseopes  as  diseussed  in  ehapter  4.  Figure  30  shows  the 
preamplifier  outputs  resulting  from  a  eoineident  event  reeorded  on  two  strips,  one  front 
strip  and  one  rear  strip.  This  event  is  the  result  of  a  gamma  depositing  energy  in  only  one 
deteetor  pixel.  It  is  elear  that  the  same  amount  of  eharge  is  eolleeted  on  eaeh  strip.  Thus 
it  is  assumed  that  the  gamma  deposited  energy  in  only  one  pixel  and  the  entire  eharge 
was  eolleeted. 
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Figure  30:  Two  strip  event  eharge  eolleetion 


Figure  31  shows  the  resulting  preamplifier  outputs  resulting  from  a  eoineident 
event  reeorded  on  three  strips,  one  rear  strip,  R3,  and  two  front  strips,  FI  and  F2. 
Although  not  obvious  at  the  seale  presented,  the  pulse  reeorded  on  strip  F2  was  reeorded 
approximately  25  ns  before  the  pulse  on  strip  FI.  This  leads  to  the  assumption  that  the 
original  gamma  entered  the  deteetor  at  the  pixel  defined  by  the  interseetion  of  strips  F2 
and  R3.  The  pulse  reeorded  on  strip  F2  was  apparently  the  result  of  one  or  more 
Compton  seattering  events  whieh  eventually  eaused  the  redueed  energy  gamma  to  enter 
the  pixel  defined  by  the  interseetion  of  strips  FI  and  R3  where  it  eontinued  to  deposit 
energy.  The  single  pulse  reeorded  by  strip  R3  is  a  result  of  the  very  short  amount  of  time 
between  the  Compton  seatter  interaotion(s)  and  the  final  photoeleetrie  absorption 
interaetion. 


60 


Volts  fmVl 
800- 

600- 


200  400  600 

Time  [/is] 


Figure  31;  Three  strip  event  eharge  eolleetion 


The  eharge  sharing  between  the  two  front  strips  is  obvious.  Close  analysis  shows 
that  the  sum  of  the  eharge  eolleeted  on  the  two  rear  strips  is  equal  to  the  eharge  eolleeted 
on  the  front  strip.  The  same  results  were  found  if  the  eharge  sharing  was  between  two 
rear  strips.  Onee  again,  this  leads  to  the  assumption  that  full  eharge  eolleetion  is  being 
aeeomplished. 


Preamplifier  Decay  Constant  (TAU) 

One  of  the  user  inputs  within  the  XIA  software  is  the  preamplifier  deeay  eonstant. 
As  diseussed  in  4,  the  DGF  software  eontains  a  funetion  for  obtaining  this  value.  It  does 
so  by  eomputing  ten  values  for  the  preamplifier  deeay  eonstant  of  a  speeified  strip  and 
returning  the  average  of  the  ten.  The  user  is  then  able  to  keep  the  eomputed  value,  obtain 
a  new  average,  or  use  a  built-in  funetion  to  find  the  value  eomputed  from  a  single  pulse. 
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For  this  work  I  utilized  the  built-in  averaging  funetion  and  obtained  the  values  shown  in 
Table  6. 


Table  6;  Pre-amplifier  decay  constants 


Strip 

Tau  [/Us] 

FI 

47.4 

F2 

43.8 

F3 

41.8 

F4 

52.1 

R1 

43.1 

R2 

45.0 

R3 

47.6 

R4 

45.7 

Energy  Calibration 

As  described  in  Chapter  4,  proper  energy  calibration  is  crucial  in  order  to  image 
within  a  given  energy  range.  This  requires  that  the  gain  setting  for  each  strip  be  adjusted 
so  that  the  centroid  of  a  specified  full  energy  peak  is  recorded  in  a  predetermined 
channel.  For  this  work,  the  specified  peak  centroid  of  392  keV  was  aligned  with  channel 
20200.  As  shown  in  Table  7,  this  required  different  gain  settings  for  each  strip.  Because 
of  the  inability  to  reproduce  exact  high  voltage  bias  values,  new  gain  setting  adjustments 
are  required  when  the  high  voltage  bias  is  removed  from  the  detector.  The  settings 
shown  below  are  those  used  when  determining  the  energy  resolution  of  the  detector. 

Each  setting  was  verified  by  conducting  several  runs  to  ensure  the  centroid  channel 
number  did  not  shift  by  more  than  +  1  channel.  According  to  these  settings,  the  largest 
difference  between  a  particular  centroid  on  different  charge  collection  strips  should  be 
four  channels. 
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Table?:  Gain  settings 


Output 

Gain 

Channel 

El 

1.3084 

20203 

E2 

1.2718 

20202 

E3 

1.1970 

20198 

E4 

1.1967 

20200 

R1 

1.2895 

20202 

R2 

1.2347 

20201 

R3 

1.2460 

20199 

R4 

1.2225 

20200 

After  all  eight  eharge  eolleetion  strips  had  been  ealibrated,  the  energy  ealibration 
had  to  be  known  for  use  in  the  image  proeessing  eode.  The  energy  of  an  event  in  raw 
data  is  reeorded  in  the  list  mode  data  eolleetion  as  a  ehannel  number.  There  is  no  way  of 
determining  the  energy  ealibration  seale  within  the  XIA  software.  So  the  raw  data  were 
read  into  the  image  proeessing  eode  as  ehannel  numbers,  and  the  eode  allows  the  user  to 
either  aeeept  a  default  ealibration  or  manually  enter  the  energy  ealibration.  The  default 
energy  ealibration  used  by  the  image  proeessing  eode  used  six  full  energy  peaks,  whose 
values  of  ehannel  and  energy  were  plotted  in  Exeel.  Onee  the  data  were  plotted.  Exeel 
performed  a  least  squares  fit  on  the  data  and  determined  an  equation  of  the  line.  The  plot 
of  the  data  is  shown  in  Eigure  32. 
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Figure  32:  Energy  ealibration  of  the  HPGeDSSD 


The  ealeulated  value  of  the  least  squares  fit  is 

Channel  =  38.852  *  Energy[keV]  +  4979.8  (6) 

Sinee  the  data  were  eolleeted  in  MCA  mode  with  an  energy  eutoff  of  5000,  the  data  were 
shifted  to  the  right  by  5000  ehannels.  This  does  not  oeeur  when  operating  in  list  mode 
data  eolleetion,  so  the  energy  fit  used  within  the  eode  is 

Channel  =  38.852  *  Energy[keV]  -  20.2  (7) 


Rise  Time  Optimization 

The  resolution  of  the  deteetor  is  dependent  on  both  the  rise  time  and  flat  top  time 
of  the  energy  trapezoidal  filter.  The  XIA  doeumentation  suggests  that  the  resolution  is 
relatively  independent  of  the  flat  top  time,  and  highly  dependent  on  the  energy  filter  rise 
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time.  It  also  suggests  that  the  values  of  the  energy  filter  rise  time  should  be  varied  in 
order  to  determine  the  optimal  rise  time  for  each  charge  collection  strip.  To  optimize  the 
rise  time  of  individual  strips,  eleven  values  of  the  energy  fdter  rise  time  were  used,  and 
the  FWHM  of  the  392  keV  full  energy  peak  of  a  tin  strip  source  were  measured.  Figure 
33  is  a  plot  of  the  eleven  values  obtained  for  each  strip. 
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Figure  33;  Plot  of  FWFIM  vs.  energy  filter  rise  time 


The  minimum  resolution  for  each  strip  occurred  at  various  rise  time  values. 
Therefore,  the  best  optimization  should  be  achievable  by  using  the  best  rise  time  for  each 
charge  collection  strip  individually.  However,  because  of  a  limitation  in  the  DGF-4C 
controlling  software,  using  various  energy  filter  rise  times  restricts  the  coincident  time 
window  allowed  for  determining  coincident  events  in  adjacent  strips.  The  desired 
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coincidence  window  is  100  ns.  If  two  strips  have  energy  filter  rise  times  varying  by  just 
0.4  yus,  the  shortest  coincidence  window  allowed  is  400  ns.  This  led  to  the  use  of  a  2.4  //s 
energy  filter  rise  time  for  all  charge  eollection  strips.  This  value  was  ehosen  in  order  to 
maintain  the  narrowest  energy  window  for  imaging. 


Energy  Resolution  of  the  Detector 

In  order  to  eharaeterize  a  deteetor,  the  energy  resolution  must  be  known.  Chapter 
4  discusses  the  process  of  measuring  the  energy  resolution.  Since  the  HPGeDSSD  has 
ten  eharge  readouts,  each  of  these  must  be  charaeterized.  Figure  34  is  one  example  of  the 
spectra  collected  to  measure  the  resolution.  All  eight  charge  collection  strips  were 
measured  and  the  remaining  seven  spectra  are  shown  in  Appendix  C.  The  time  for  the 
data  collection  was  36000  seconds. 


Figure  34;  Strip  R4  activated  Sn  spectrum  used  for  resolution  measurement 

From  eaeh  speetrum  the  XIA  software  was  used  to  find  the  absolute  FWHM  of 
each  of  the  full  energy  peaks.  Figure  35  shows  the  FWHM  measurements  for  all  of  the 
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charge  collection  strips  using  a  2.4  //s  energy  filter  rise  time.  These  measurements  of 
FWHM  can  be  used  to  help  determine  the  width  of  the  energy  window  required  for 
imaging.  The  energy  window  for  imaging  a  particular  isotope  must  be  wide  enough  to 
include  the  minimum  and  maximum  expected  values  for  the  particular  full  energy  peak 
across  all  charge  collection  strips.  For  this  work  it  was  found  that  an  energy  window  of 
about  6  keV  is  sufficient  for  comparing  charge  pulses  collected  on  opposite  faces  of  the 
detector.  This  value  corresponds  to  approximately  two  times  the  FWFIM  for  the  worst 
charge  collection  strip,  namely  R2. 
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Figure  35:  Plot  of  FWHM  vs.  energy  for  individual  charge  collection  strips 
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Peak  Efficiency  of  the  Detector 


The  peak  effieieney  of  the  deteetor  must  be  known  in  order  to  determine  the  time 
required  to  produee  an  image.  As  diseussed  in  Chapter  4,  the  peak  effieieney  was 
measured  for  six  full  energy  peaks.  Table  8  shows  the  values  of  effieieney  that  were 
measured  using  four  radiation  sourees,  and  Figure  36  shows  the  data  plotted  on  a  log-log 
seale. 


Table  8:  Deteetor  effieieney  measurements 


Full 

energy 

peak 

[keV] 

Intrinsie 

Peak 

Effieieney 

122 

0.44121 

244 

0.088358 

344 

0.042821 

392 

0.031591 

511 

0.020091 

662 

0.010875 

Energy  [keV] 


Figure  36:  Deteetor  peak  effieieney  measurements 
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The  higher  effleieney  found  around  122  keV  is  due  to  the  faet  that  photoeleetrie 
absorption  dominates  over  Compton  seattering  up  to  about  200  keV.  This  domination  of 
photoeleetrie  absorption  inereases  the  likelihood  that  these  lower  energy  gammas  will 
deposit  all  of  their  energy  within  the  germanium.  On  the  other  hand,  as  the  energy 
inereases  the  likelihood  inereases  that  the  gamma  will  undergo  one  or  multiple  Compton 
seatters  and  then  depart  from  the  deteetor.  This  results  in  a  lower  pereentage  of  gammas 
that  deposit  their  entire  energy  in  the  germanium,  thus  deereasing  effleieney  with 
inereasing  gamma  energy.  More  analysis  would  be  required  to  determine  an  optimum 
deteetor  erystal  design  whieh  would  inerease  the  peak  effleieney  of  the  deteetor. 


Investigation  of  edge  pixel  collimator  holes 

As  discussed  in  Chapter  4, 1  investigated  the  feasibility  of  using  collimator  holes  located 
at  the  edge  of  the  detector  pixels.  Other  research  investigating  the  response  near  a  strip 
gap  was  performed  at  the  European  Synchrotron  Radiation  Facility  [12].  This  work  used 
a  source  beam  of  x-rays  ranging  from  15  to  100  keV  collimated  to  10  microns  wide.  The 
results  showed  that  a  very  large  number  of  full  energy  peaks  were  recorded  when  the 
source  was  about  50  microns  from  the  gap.  This  would  then  require  a  collimator  hole 
width  of  less  than  100  microns  to  ensure  that  charge  would  be  shared  between  two  strips 
most  of  the  time.  What  1  had  hoped  to  And  was  that  the  distance  from  the  gap  which 
resulted  in  a  large  percentage  of  three  strip  events  would  be  signiflcantly  larger  for  higher 
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energy  photons,  resulting  in  the  ability  to  use  edge  pixel  holes  with  high  enough 
efficiency  for  imaging. 

Using  a  1  mm  by  25  mm  collimated  tin  strip  source,  I  gathered  spectrum  at 
various  locations  near  the  edge  between  strips  R2  (channel  #1)  and  R3  (channel  #  2). 
Figure  37  shows  the  spectrum  taken  with  the  center  of  the  collimated  source  located  0.5 


mm  from  the  gap  towards  the  center  of  strip  R2. 


As  the  spectrum  above  shows,  both  strips  recorded  a  very  large  number  of  full 
energy  peaks.  The  current  system  of  processing  an  event  does  not  provide  a  way  to 
distinguish  a  full  energy  peak  recorded  near  the  gap  from  a  full  energy  peak  recorded  at 
the  center  of  a  detector  strip.  This  causes  all  full  energy  peaks  to  be  recorded  as  center 
pixel  events.  Only  full  energy  peaks  which  share  charge  between  two  strips  on  one  side 
of  the  detector  and  is  recorded  as  a  full  energy  peak  on  the  opposite  side  could  be 
assigned  to  an  edge  pixel.  Therefore,  the  edge  pixel  hole  diameter  would  have  to  be  less 
than  2  mm.  A  hole  less  than  2  mm  diameter  does  not  provide  enough  efficiency  for 
imaging.  Also,  in  order  to  maintain  a  true  image,  all  holes  should  be  the  same  size  in 
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order  to  provide  approximately  the  same  effieieney.  This  would  preelude  the  use  of 
different  hole  sizes  for  edge  pixels  and  eenter  pixels,  whieh  would  require  obtaining  a 
good  point  spread  funetion  for  eaeh  pixel  based  on  its  size  and  loeation. 

The  results  of  the  edge  effeet  test  showed  that  I  would  not  able  to  inerease  the 
spatial  resolution  of  the  system  by  adding  edge  hole  pixels.  This  led  me  to  use  only  one 
eollimator  hole  per  deteetor  pixel.  However,  I  was  still  able  to  use  three  strip  events  by 
assigning  them  to  a  single  deteetor  pixel. 


Collimator  Alignment 

Proper  alignment  of  the  eonverging  eollimator  is  eritieal  for  proper  spatial 
resolution  of  an  image.  As  diseussed  in  Chapter  4,  the  first  step  in  aligning  the  eollimator 
was  to  align  the  eenter  hole  of  the  eollimator  with  the  eenter  pixel  of  the  deteetor.  This 
was  done  by  aligning  the  eollimator  with  the  approximate  eenter  of  the  deteetor  using  a 
Eu-152  souree  on  the  eenter  eollimator  hole.  Speetra  were  then  taken  and  the  eollimator 
position  was  adjusted  aeeording  to  the  speetra.  This  proeess  was  repeated  until  the 
souree  was  seen  only  at  the  eenter  pixel.  Figure  38  and  Figure  39  show  the  final  speetra 
taken  for  this  step.  Sinee  the  peaks  were  elearly  present  only  on  strips  R3  and  F3,  both 
ehannel  #2,  it  was  assumed  that  the  eollimator  was  indeed  eentered  on  the  deteetor. 
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The  next  step  in  aligning  the  eollimator  involved  plaeing  the  source  at  the  top 
right  image  pixel  of  the  detector  and  taking  spectra.  Just  as  before,  the  spectra  were 
analyzed  and  the  collimator  was  adjusted  accordingly.  After  only  a  few  adjustments,  the 
top  right  pixel  was  properly  aligned.  Figure  40  and  Figure  41  show  the  final  spectra 
taken,  which  verified  proper  alignment  of  the  top  right  pixel.  Proper  alignment  is 
inferred  from  the  fact  that  on  strips  R1  and  FI  show  the  peak. 
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Counts 


The  same  proeess  was  repeated  for  the  bottom  left  image  pixel  to  verify  proper 
alignment.  Figure  42  and  Figure  43  show  the  resulting  speetra  for  the  bottom  left  pixel. 
Just  as  before,  proper  alignment  was  inferred  from  the  faet  that  strips  R4  and  F4  show  the 
peak.  With  the  deteetor  properly  aligned,  I  was  then  ready  to  begin  taking  image  data. 
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Figure  42;  Rear  strip  spectra  for  bottom  left  pixel  alignment 
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Figure  43;  Front  strip  spectra  for  bottom  left  pixel  alignment 
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Image  Formation 


Several  representative  images  were  eollected  to  demonstrate  the  performance  of 
the  detector.  As  discussed  in  Chapter  4,  the  areas  of  interest  for  image  analysis  included 
the  utilization  of  two  strip  events  versus  three  strip  events,  the  ability  to  distinguish  a 
point  source  from  a  distributed  source,  and  the  ability  to  image  using  multiple  energy 
windows.  The  resulting  MATLAB  outputs  are  shown  below.  Each  output  consists  of  a 
four  by  four  pixel  standard  grayscale  image,  a  three  dimensional  surface  plot  of  the  image 
showing  relative  counts  in  each  pixel,  and  a  histogram  showing  the  number  of  counts  for 
each  1  keV  increment  within  the  specified  energy  windows.  For  each  of  the  standard 
pixel  images  and  surface  plots,  black  represents  zero  counts  and  white  indicates  the 
highest  number  of  counts  within  the  image.  A  grayscale  legend  is  also  provided  to  assist 
in  comparing  variations  in  the  image. 

Figure  44  shows  the  image  results  of  a  Csl37  point  source  that  was  imaged  for 
9426  seconds  and  processed  using  only  two  strip  events  and  a  660-666  keV  energy 
window.  Table  9  shows  the  resulting  image  matrix  for  Figure  44. 
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Surface  plot  of  660-666  keV  660-666  keV  image 


Energy  [keV] 

Figure  44;  Cs-137  point  source  image  results  using  only  two  strip  events 


Table  9;  Cs-137  point  source  image  matrix  using  only  two  strip  events 
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Good  collimation  of  the  source  in  Figure  44  is  evident  by  the  relatively  low 
number  of  counts  recorded  in  pixels  which  were  not  visible  to  the  source.  Several  other 
pixels  were  tested  for  good  collimation  with  similar  results  to  those  shown  above.  Time 
precluded  gathering  enough  data  to  develop  a  point  spread  function  to  account  for 
efficiency  differences  between  the  detector  pixels.  Therefore,  for  this  work  it  was 
assumed  that  the  efficiency  from  pixel  to  pixel  was  relatively  constant. 
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Figure  45  shows  the  image  results  of  processing  the  same  data  set  as  that  used  in 
Figure  44  using  both  two  and  three  strip  events  and  the  same  660-666  keV  energy 
window.  Table  10  shows  the  resulting  image  matrix  for  Figure  45. 


Surface  plot  of  660-666  keV 
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Figure  45:  Cs-137  point  source  image  results  using  two  and  three  strip  events 


Table  10:  Cs-137  point  source  image  matrix  using  two  and  three  strip  events 
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Comparison  of  Figure  44  and  Figure  45  shows  that  the  addition  of  three  strip 
events  increased  number  of  counts  recorded  in  the  illuminated  pixel  by  28%.  On  the 
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other  hand,  this  increase  in  counts  comes  at  the  cost  of  significantly  increasing  the 
number  of  counts  in  the  adjacent  nearest  neighbor  pixels.  At  first  glance  one  may 
conclude  that  the  addition  of  three  events  degrades  the  image.  However,  many  more 
images  would  be  required  in  order  to  quantify  the  effect  of  adding  three  strip  events.  For 
applications  which  involve  very  low  count  rate,  such  as  with  on  site  weapon  inspections, 
the  addition  of  three  strip  events  could  prove  crucial  to  obtaining  enough  counts  for  a 
viable  image.  Further  research  is  required  to  determine  when  the  use  of  three  strip  events 
can  be  justified. 

Figure  46  shows  the  image  results  of  an  activated  Sn  distributed  source  comprised 
of  five  Sn  strips  in  the  shape  of  the  letter  “A”  which  was  imaged  for  30526  seconds, 
approximately  8.5  hours.  The  image  was  developed  using  two  energy  windows,  namely 
389-395  keV  and  424-430  keV.  Recall  that  these  are  the  peaks  chosen  to  simulate  Pu- 
239  in  a  weapon  storage  cask.  The  histogram  of  the  image  clearly  shows  that  each 
energy  window  made  a  significant  contribution  to  the  final  image.  Figure  47  and  Figure 
48  show  the  contributions  each  energy  window  made  in  Figure  46.  These  images  show 
the  ability  to  image  a  distributed  source  of  SNM  using  either  a  single  or  multiple  energy 
windows.  Comparison  of  these  distributed  source  images  to  those  obtained  from  a  point 
source  clearly  show  the  ability  of  the  system  to  distinguish  a  point  source  from  a 
distributed  source. 
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389-395  and  424430  keV  image 
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Figure  46:  Image  results  of  distributed  Sn  souree  using  two  energy  windows 
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Figure  47:  Image  results  of  distributed  Sn  souree  using  a  389-395  keV  energy  window 
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424430  keV  image 
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Figure  48;  Image  results  of  distributed  Sn  souree  using  a  424-430  keV  energy  window 


Close  analysis  of  the  three  figures  above  shows  that  the  best  qualitative  image 
eontrast  was  aetually  achieved  using  only  the  389-395  keV  energy  window.  This  may 
lead  one  to  assume  that  you  really  only  need  to  image  in  one  energy  window  at  a  time. 
However,  keep  in  mind  that  one  major  advantage  in  being  able  to  image  in  multiple 
energy  windows  comes  from  the  fact  that  weapon  inspection  times  are  limited,  and  that 
the  flux  of  measurable  gamma  photons  emitted  through  a  weapon  cask  is  quite  small. 
Therefore,  having  the  ability  to  combine  counts  from  numerous  peaks  greatly  increases 
the  overall  efficiency  of  the  system.  The  image  processing  code  could  easily  be  updated 
to  handle  as  many  energy  windows  as  desired. 

Another  way  in  which  imaging  across  multiple  energy  windows  can  increase  the 
efficiency  of  the  system  is  by  providing  the  ability  to  image  multiple  nuclides 
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simultaneously  or  individually  using  a  single  data  set.  The  presence  and  physical 
orientation  of  multiple  isotopes  could  provide  information  about  the  weapon  design. 
However,  some  caution  must  be  used  when  using  multiple  energy  windows  to  image 
multiple  nuclides.  The  key  lies  in  the  fact  that  the  highest  value  recorded  on  the  image  is 
the  highest  count  for  any  pixel  within  the  specified  energy  windows.  An  example  of  this 
is  shown  in  Figure  49.  The  image  is  of  a  Eu-152  point  source  in  the  bottom  left  comer 
and  a  Cs-137  point  source  near  the  top  right  comer  imaged  across  the  341-347  keV  (Eu- 
152)  and  660-665  keV  (Cs-137)  energy  windows. 


Surface  plot  of  341-347  and  660-665  keV  341-347  and  660-665  keV  image 


Energy  [keV] 

Eigure  49;  Two  source  image  of  Eu-152  and  Cs-137  point  sources 

The  higher  activity  of  the  Eu-152  source  is  evident  by  the  higher  number  of 
counts  recorded  under  the  341-347  keV  window.  Careful  choices  of  the  energy  windows 
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clearly  make  visible  the  orientation  of  the  two  sources.  However,  it  is  not  elear  if  the 
visible  shapes  consist  of  a  mixture  of  the  two  nuclides  or  if  they  are  in  fact  two  separate 
sources.  Additional  images  ean  be  produced  using  the  same  data  set  with  individual 
energy  windows  to  ensure  that  in  fact  the  two  apparent  sources  are  different  nuclides  and 
that  the  eounts  for  a  partieular  pixel  is  not  from  a  mixture  of  the  two  nuclides.  Figure  50 
and  Figure  51  show  the  resulting  images  using  the  energy  windows  for  each  source 
individually.  The  resulting  image  in  Figure  50  is  aetually  corrupted  by  higher  energy  Cs- 
137  gammas  whieh  are  Compton  scattered  down  into  the  energy  window  being  viewed. 


Surface  plot  of  341-347  keV  341-347  keV  image 
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Figure  50:  Two  source  image  using  341-347  keV  (Eu-152)  energy  window 
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Surface  plot  of  660-665  keV  660-665  keV  image 
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Figure  51:  Two  source  image  using  660-665  keV  (Cs-137)  energy  window 

Several  other  multi-nuclide  images  were  developed  during  this  work  with  results 
comparable  to  those  presented  above.  This  analysis  resulted  in  a  high  degree  of 
confidence  of  the  use  of  multiple  energy  windows  to  identify  multiple  nuclides  within  a 
source  image. 

The  ability  of  the  system  to  distinguish  a  distributed  source  from  a  point  source  in 
a  configuration  consistent  with  that  encountered  during  on  site  weapons  inspections  was 
tested  using  the  pit  simulation  source  described  in  Chapter  3  and  a  Eu-152  point  source. 
Each  source  was  placed  71  cm  from  the  collimator  along  the  center  axis  of  the  4  by  4 
pixel  array  used  for  imaging.  Eigure  52  shows  the  resulting  image  of  the  pit  simulation 
source  for  a  collection  time  of  14964  seconds.  Although  a  circular  symmetric  image  was 
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not  obtained,  it  is  still  clear  that  the  source  is  indeed  distributed.  The  non-symmetric 
image  could  be  a  result  of  the  source  leakage  discussed  in  Chapter  3. 


Surface  plot  of  660-667  keV  660-667  keV  image 
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Figure  52:  Image  of  Cs-137  pit  simulation  source  71  cm  from  detector 


Figure  53  shows  the  image  acquired  from  the  Eu-152  source  at  a  distance  of  71 
cm  for  a  collection  time  of  approximately  10  hours.  Comparison  of  Figure  52  and  Figure 
53  clearly  shows  the  ability  of  the  system  to  distinguish  a  point  source  from  a  distributed 
source  in  a  configuration  consistent  with  that  encountered  during  on  site  weapon 
inspections. 
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Surface  plot  of  341-347  keV 
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Figure  53:  Image  of  Eu-152  souree  71  em  from  deteetor 
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VI,  Conclusions  and  Recommendations 


Figure  54;  Combination  of  four  images  to  spell  DTRA 


Figure  54  was  ereated  by  eombining  four  images  that  were  produeed  by  plaeing 
the  aetivated  Sn  strip  sourees  in  a  eonfiguration  that  formed  the  letters  D,T,R,  and  A.  An 
example  of  the  placement  of  the  radiation  sources  for  the  “T”  is  shown  in  Figure  55. 


Chapter  Overview 

The  principle  objective  of  this  project  was  to  develop  an  improved  FlPGe  strip 
detector  imaging  system  to  help  determine  the  basic  size,  shape  and  composition  of  SNM 
in  a  nuclear  weapon  pit.  Key  areas  of  interest  included  the  performance  of  the  detector, 
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performance  of  the  input  logic  module,  performance  of  the  digital  gamma  finder 
modules,  the  application  of  a  converging  hole  collimator,  and  application  of  a  new  image 
processing  code. 

Detector  Performance 

The  performance  of  the  detector  was  essentially  what  was  expected  for  a  high 
purity  germanium  detector.  The  resolution  of  the  charge  collecting  strips  was  slightly 
poorer  than  a  typical  HPGe  detector,  but  the  efficiency  measurements  were  in  agreement 
with  expectations  for  a  detector  of  this  size.  The  XIA  module  control  software  version 
3.04  proved  to  be  an  improvement  over  version  3.02  which  was  used  in  previous  work 
with  this  detector. 

Resolution;  The  resolution  measurements  made  using  the  XIA  modules  were 
consistently  poorer  (larger)  than  those  quoted  by  ORTEC,  the  detector  manufacturer. 
Table  1 1  shows  the  resolutions  reported  by  the  manufacturer,  the  measured  resolutions 
using  the  XIA  hardware  with  optimized  and  non-optimized  settings,  and  the  measured 
resolutions  using  the  settings  required  for  imaging.  As  mentioned  in  Chapter  4,  the 
energy  filter  rise  time  had  to  be  the  same  value  for  all  charge  collection  strips  in  order  to 
maintain  a  100  ns  coincidence  window  for  imaging.  This  explains  why  the  imaging 
settings  provided  optimum  resolution  for  only  two  strips,  R1  and  R2.  Further 
investigation  is  needed  to  determine  if  different  module  settings  can  be  used  to  truly 
optimize  the  resolution  of  all  detector  strips  at  the  same  time  while  not  sacrificing  the 
length  of  the  coincidence  window. 
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Table  1 1 ;  Detector  FWHM  measurements 


Strip 

Manufacturer 
Reported  FWHM 
(  122  keV) 

XIA 

Optimized 
Settings 
(635  keV) 

XIA 

Non-optimized 
Settings 
(635  keV) 

Imaging  Settings 
(635  keV) 

FI 

1.64 

2.53 

3.08 

2.76 

F2 

1.62 

2.47 

3.07 

2.78 

F3 

1.23 

2.11 

2.80 

2.55 

F4 

1.31 

2.31 

2.53 

2.41 

R1 

1.48 

2.44 

2.61 

2.44 

R2 

1.50 

2.81 

3.37 

2.81 

R3 

1.41 

2.25 

2.65 

2.56 

R4 

1.40 

2.24 

2.79 

2.76 

Efficiency:  The  NRL  has  a  45mm  x  45mm  x  10mm  germanium  strip  detector, 
with  an  efficiency  of  2%  at  662  keV.  The  detector  efficiency  calculated  during  this  work 
for  a  25mm  x  25mm  x  9mm  germanium  strip  detector  was  1.0875%  at  662  keV. 
Considering  the  difference  in  detector  thickness  as  viewed  from  the  source,  the  results 
from  this  work  agree  well.  Also  the  plot  of  efficiency  as  a  function  of  energy,  shown  in 
Figure  36,  followed  the  expected  trend  of  decreasing  efficiency  with  increasing  gamma 
energy. 


Input  Logic  Module  Performance 

The  inability  to  start  and  stop  DGF  modules  at  the  same  time  during  previous 
research  with  this  system  was  overcome  by  the  addition  of  the  ORTEC  input  logic 
module.  The  input  logic  module  provided  the  appropriate  signal  to  start  and  stop  both 
modules  at  the  same  time.  This  was  evident  by  the  ability  to  accurately  pair  together 
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events  eoineident  in  both  time  and  energy  whieh  oecurred  after  both  modules  were 
stopped  and  restarted  by  the  input  logic  module. 

Digital  Gamma  Finder  Module  Performance 

The  DGF  modules  proved  to  be  quite  effective  for  this  work.  The  ability  to 
change  settings  for  individual  strips  and  individual  modules  was  crucial  to  optimizing  the 
system.  Being  able  to  analyze  data  using  the  DGF  controlling  software  greatly  reduced 
the  amount  of  time  required  to  analyze  data  and  make  any  necessary  adjustments.  The 
biggest  shortfall  of  the  modules  is  the  inability  to  download  information  in  real  time. 

Collimator  Performance 

The  converging  hole  collimator  designed  and  built  for  this  research  effectively 
collimated  all  sources  used  for  imaging.  While  imaging  point  sources,  the  resulting 
images  had  relatively  very  few  counts  recorded  in  outlying  pixels  which  were  not 
illuminated  by  the  sources.  As  expected,  while  imaging  the  distributed  sources,  more 
events  were  recorded  in  the  outlying  pixels  than  with  the  point  sources.  However,  the 
collimation  was  sufficient  to  produce  images  with  good  contrast  between  illuminated  and 
non-illuminated  pixels. 

The  collimator  was  designed  to  provide  the  ability  to  image  sources  up  to  a  5  cm 
square  shape  at  a  distance  of  7 1  cm.  This  ability  was  demonstrated  by  the  use  of  the  pit 
simulation  source. 
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Imaging  Performance 


The  images  produeed  during  this  work  met  all  expeetations  of  being  able  to 
distinguish  a  point  souree  from  a  distributed  souree.  The  eombination  of  the  new  image 
proeessing  eode  and  the  eonverging  eollimator  made  it  possible  to  distinguish  a  point 
souree  from  a  distributed  souree  even  in  a  eonfiguration  eonsistent  with  those 
eneountered  during  on  site  weapon  inspeetions.  The  new  image  proeessing  eode  also 
provided  the  ability  to  image  aeross  multiple  energy  windows.  This  allows  the  user  to 
produee  a  viable  image  in  shorter  time  by  adding  together  the  eounts  within  multiple 
energy  windows  for  a  single  souree  nuelide.  It  also  provides  the  ability  to  obtain 
information  about  the  relative  position  of  multiple  nuelides  within  the  souree  region. 

An  additional  improvement  of  the  image  proeessing  eode  was  the  deerease  in  the 
amount  of  time  to  proeess  a  single  image.  The  previous  eode  required  several  hours  to 
produee  an  image  aeross  a  single  energy  window.  The  new  eode  ean  proeess  an  image 
aeross  two  energy  windows  in  approximately  30  seeonds.  This  greatly  reduees  the 
amount  of  time  required  for  image  analysis. 


Recommendations  for  Future  Work 

The  eurrent  XIA  software  does  not  allow  data  to  be  downloaded  to  the  host 
eomputer  until  the  entire  data  set  has  been  eolleeted.  It  would  have  been  extremely 
helpful  if  the  data  eould  have  been  transferred  to  the  host  eomputer  at  near  real  time 
This  would  have  allowed  an  image  to  be  built  up  over  time  and  the  operator  eould 
determine  when  enough  data  had  been  eolleeted.  Currently,  an  estimate  of  the  time 


91 


required  must  be  made  and  the  data  eolleeted  followed  by  image  proeessing  to  determine 
if  there  was  indeed  enough  time. 

Beeause  of  limitations  on  the  number  of  inputs  for  eaeh  DGF  module,  only  a  4  by 
4  pixel  image  was  produeed.  Although  this  was  suffieient  to  distinguish  a  point  souree 
from  a  distributed  souree,  the  addition  of  a  third  DGF  module,  along  with  some  basie 
updates  to  the  image  proeessing  eode,  would  make  it  possible  to  produee  a  5  by  5  pixel 
image  with  the  eollimator  built  for  this  work.  More  image  pixels  over  the  same  area 
would  theoretieally  allow  one  to  obtain  better  shape  information. 

A  more  important  eontribution  whieh  ean  be  made  to  this  projeet  is  developing 
the  ability  to  get  sub  pixel  resolution  for  the  eharge  eolleeting  strips.  The  idea  of  using 
edge  pixel  eollimator  holes  proved  to  be  ineffieient  at  high  photon  energies.  However, 
aeeording  to  the  Navy  Researeh  Lab  [8],  it  is  possible  to  determine  whieh  portion  of  the 
eharge  eolleeting  strip  aetually  eolleeted  the  eharge  by  viewing  the  indueed  noise  of  the 
two  strips  adjaeent  to  the  one  that  eolleeted  the  eharge.  This  would  allow  an  image  array 
of  10  X  10,  if  the  position  eould  be  determined  within  2.5  mm.  The  advantage  to 
proeessing  the  data  in  this  way  to  inerease  the  resolution  as  opposed  to  simply  ereating  a 
strip  deteetor  with  10  strips  on  eaeh  side  is  that  the  inerease  in  hardware  is  minimal  with 
this  method,  and  if  the  system  is  to  be  used  for  treaty  monitoring,  portability  is  a  eoneern 
and  hardware  must  be  minimized. 

Summary 

The  ability  to  use  an  HPGeDSSD  eoupled  with  a  eonverging  hole  eollimator  as  an 
imaging  deviee  to  deteet  spoof  weapons  has  been  demonstrated.  Although  the  spatial 
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resolution  of  the  deteetor  eould  not  be  improved  using  edge  hole  pixels,  the  eurrent 
spatial  resolution  proved  to  be  suffieient  for  distinguishing  a  point  souree  from  a 
distributed  souree.  The  use  of  a  eonverging  eollimator  provided  the  ability  to  image  a 
souree  that  was  larger  than  the  deteetor  in  a  eonfiguration  eonsistent  with  that 
eneountered  during  on  site  weapon  inspeetions.  The  addition  of  an  input  logie  module, 
eombined  with  a  new  image  proeessing  eode,  made  it  possible  to  aequire  and  proeess 
three  strip  events  for  imaging.  The  new  image  proeessing  eode  also  provided  the  ability 
to  simultaneously  identify  multiple  nuelides  by  imaging  aeross  multiple  energy  windows. 
In  addition,  this  was  all  made  possible  while  greatly  redueing  the  time  required  for 
proeessing  an  image. 
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Appendix  A.  Collimator  Design  Code 


Graphics'Graphics' 

Collimator  design  code 

Purpose:  This  code  is  used  to  help  design  efective  collimators  for  imaging  special 
nuclear  material  using  a  HPGe  Double  Sided  Strip  Detector  with  a  5  x  5  array  of  5  x  25 
mm  strips. 

CPT  Doug  Rothenbush 
Version  1.0 
6Dec04 

All  measurements  are  in  [mm]  unless  noted. 

Input  parameters 

holewidth=2.9; 

collimatorthickness=190; 

focallength=1160; 

caskradius=710; 

Find  slope  of  centerlines 


lineslcpe[l]  = 


50. 

focallength 


angle[l]  =  Ar«a:an[lineslcpe[l]  ] 


180 

*  -  r 

TT 


lineslcpe[2]  = 


25. 

focallength 


angle[2]  =  Atxiran[lineslcpe[2]  ] 


180 

*  -  r 

TT 


lineslcpe[3]  =  0 


angle[3]  =  0; 

0.0431034 

0.0215517 

0 

Determine  critical  locations 

sourceface=focallength-caskradius  (*  collimator  must  be  atleast  the  cask  radius  from  the 
source  *) 

450 

detface=sourceface-collimatorthickness 

260 

10 

crystalcenter  =  - 

lineslqpe[l] 


232. 


detectorhousmg=crystalcenter+19.5  (*  must  be  less  than  detface  *) 
251.5 


Compute  the  y  value  of  the  centerlines  at  the  collimator  source  face  (inside  face  of  source 
endplate)  [mm] 

Do[sourcecenter[n]=sourceface*lmeslope[n],{n,l,3}] 

sourcecenter[l] 

sourcecenter[2] 

sourcecenter[3] 
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19.3966 

9.69828 

0 

Convert  to  inches 

Do[sourcecenterinches[n]=sourcecenter[n]/25 .4,  {n,  1 ,3  }  ] 
sourcecenterinches  [  1  ] 
sourcecenterinches  [2] 
sourcecenterinches  [3  ] 

0.763644 

0.381822 

0 

Compute  the  y  value  of  the  centerlines  at  the  outside  face  of  source  endplate[mm] 

Do[outsidesourcecenter[n]=(sourceface+6.35)*lineslope[n],{n,l,3}] 

outsidesourcecenter[  1  ] 
outsidesourcecenter[2] 
outsidesourcecenter[3] 

19.6703 

9.83513 

0 

Convert  to  inches 

Do[outsidesourcecenterinches[n]=outsidesourcecenter[n]/25 .4,  {n,  1 ,3  }  ] 

outsidesourcecenterinches[  1  ] 

outsidesourcecenterinches[2] 

outsidesourcecenterinches[3] 

0.77442 

0.38721 

0 

Compute  the  y  value  of  the  centerlines  at  the  collimator  detector  face  (inside  face  of 
detector  endplate)  [mm] 

Do[detcenter[n]=detface*lineslope[n],  {n,  1 ,3 }  ] 

detcenter[l] 

detcenter[2] 

detcenter[3] 

11.2069 

5.60345 

0 

Convert  to  inches 

Do[detcenterinches[n]=detcenter[n]/25 .4,  (n,  1 ,3  }  ] 

detcenterinches[  1  ] 

detcenterinches[2] 

detcenterinches[3] 

0.441216 

0.220608 

0 

Compute  the  y  value  of  the  centerlines  at  the  outside  face  of  detector  endplate[mm] 
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Do[outsidedetcenter[n]=(detface-6.35)*lineslope[n],{n,l,3}] 

outsidedetcenter[  1  ] 
outsidedetcenter[2] 
outsidedetcenter[3] 

10.9332 

5.46659 

0 

Convert  to  inches 

Do[outsidedetcenterinches[n]=outsidedetcenter[n]/25 .4,  {n,  1 ,3 }  ] 
outsidedetcenterinches  [  1  ] 
outsidedetcenterinches  [2] 
outsidedetcenterinches  [3  ] 

0.430441 

0.21522 

0 

Find  the  width  of  the  holes  at  the  collimator  faces  and  endplate  faces 

,  .  5  *  holewidth 

Ay[l]  =  - 

Cos  [  ArcTan  [  lineslc^  [  1]  ]  ] 

.  5  *  holewidth 

Ay[2]  =  - 

Cos [ ArcTan [lineslcpe[ 2]  ]  ] 

^y[3]  =  .5*holewidth 

1.45135 

1.45034 

I. 45 

Convert  to  inches 

yinches[l]  =  y[l]/25.4 
yinches[2]  =  y[2]/25.4 
yinches[3]  =  y[3]/25.4 
0.0571396 
0.0570999 
0.0570866 

Find  y  value  for  top  of  holes  on  the  source  side  (inside  of  source  endplate) 

Do[topsource[n]=sourcecenter[n]+  y  [n] ,  {n,  1 ,3  }  ]; 

topsource[l] 

topsource[2] 

topsource[3] 

20.8479 

II. 1486 
1.45 

Convert  to  inches 

Do[topsourceinches[n]=topsource[n]/25 .4,  {n,  1 ,3  }  ] ; 

topsourceinches[  1  ] 

topsourceinches[2] 

topsourceinches[3] 

0.820783 
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0.438922 

0.0570866 

Find  y  value  for  bottom  of  holes  on  the  source  side  (inside  of  source  endplate) 
Do[bottomsource[n]=sourcecenter[n]-  y [n] ,  {n,  1 ,3 } ]; 
bottomsource[l] 
bottomsource[2] 
bottomsource[3] 

17.9452 

8.24794 

-1.45 

Convert  to  inches 

Do[bottomsourceinches[n]=bottomsource[n]/25.4,{n,l,3}]; 

bottomsourceinches[  1  ] 
bottomsourceinches[2] 
bottomsourceinches[3] 

0.706504 

0.324722 

-0.0570866 

Find  y  value  for  top  of  holes  on  the  outside  of  source  endplate 
Do[topoutsidesource[n]=outsidesourcecenter[n]-l-  y[n],{n,l,3}]; 
topoutsidesource[  1  ] 
topoutsidesource[2] 
topoutsidesource[3] 

21.1216 

11.2855 

1.45 

Convert  to  inches 

Do[topoutsidesourceinches[n]=topoutsidesource[n]/25 .4,  (n,  1 ,3  }  ] ; 

topoutsidesourceinches[  1  ] 

topoutsidesourceinches[2] 

topoutsidesourceinches[3] 

0.831559 

0.44431 

0.0570866 

Find  y  value  for  bottom  of  holes  on  the  outside  of  source  endplate 
Do[bottomoutsidesource[n]=outsidesourcecenter[n]-  y[n],{n,l,3}]; 
bottomoutsidesource[  1  ] 
bottomoutsidesource[2] 
bottomoutsidesource[3] 

18.2189 

8.38479 

-1.45 

Convert  to  inches 

Do[bottomoutsidesourceinches[n]=bottomoutsidesource[n]/25 .4,  {n,  1 ,3  }  ] ; 
bottomoutsidesourceinches[  1  ] 
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bottomoutsidesourceinches[2] 
bottomoutsidesourceinches  [3  ] 

0.71728 

0.33011 

-0.0570866 

Find  y  value  for  top  of  holes  on  the  detector  side  (inside  of  detector  endplate) 
Do[topdet[n]=detcenter[n]+  y  [n] ,  (n,  1 ,3  }  ] 
topdet[l] 
topdet[2] 
topdet[3] 

12.6582 

7.05378 

1.45 

Convert  to  inches 

Do[topdetinches[n]=topdet[n]/25 .4,  (n,  1 ,3  }  ] 

topdetinches[l] 

topdetinches[2] 

topdetinches[3] 

0.498356 

0.277708 

0.0570866 

Find  y  value  for  bottom  of  holes  on  the  detector  side  (inside  of  detector  endplate) 
Do[bottomdet[n]=detcenter[n]-  y [n],  {n,  1 ,3 } ] 
bottomdet[l] 
bottomdet[2] 
bottomdet[3] 

9.75555 

4.15311 

-1.45 

Convert  to  inches 

Do[bottomdetinches[n]=bottomdet[n]/25 .4,  {n,  1 ,3  }  ] 
bottomdetinches[  1  ] 
bottomdetinches[2] 
bottomdetinches[3] 

0.384077 

0.163508 

-0.0570866 

Find  y  value  for  top  of  holes  on  the  outside  of  detector  endplate 
Do[topoutsidedet[n]=outsidedetcenter[n]-l-  y [n] ,  {n,  1 ,3  }  ] 
topoutsidedet[l] 
topoutsidedet[2] 
topoutsidedet[3] 

12.3845 

6.91693 

1.45 
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Convert  to  inches 

Do[topoutsidedetinches[n]=topoutsidedet[n]/25 .4,  {n,  1 ,3 }  ] 
topoutsidedetinches[  1  ] 
topoutsidedetinches[2] 
topoutsidedetinches[3] 

0.48758 

0.27232 

0.0570866 

Find  y  value  for  bottom  of  holes  on  the  outside  of  detector  endplate 
Do[bottomoutsidedet[n]=outsidedetcenter[n]-  y[n], {n,  1,3}] 
bottomoutsidedet[  1  ] 
bottomoutsidedet[2] 
bottomoutsidedet[3  ] 

9.48184 

4.01626 

-1.45 

Convert  to  inches 

Do[bottomoutsidedetinches[n]=bottomoutsidedet[n]/25 .4,  {n,  1 ,3 }  ] 
bottomoutsidedetinches[  1  ] 
bottomoutsidedetinches[2] 
bottomoutsidedetinches[3] 

0.373301 

0.15812 

-0.0570866 


lin^[n,  x_] 


tx)pline[n,  x_] 


bottciiiline[n,  x  ] 


- 1  X  (x-detface)  +  tcpdet[n]  ; 


Find  the  equations  for  the  field  of  view  lines  and  the  centerlines 

r  /  bottcmsouraern]  -tx»det[n]  \ 

DoflineA[n,  x  ]  =  - )  x  (x-detface)  +  topdet[n]  ; 

'•  “  \  cx)lliiiBtQrthickiiess  ' 

I  topsource[n]  -  bottatidet[n]  \ 

lin^[n/  X  ]  =  -  X  (x-detfaoe)  +  bottaTidet[n]  ;  oaiterline[n/  x  ]  =  liJieslope[n]  *x; 

“  \  cx>lliiiBtorthickiiess  '  ~ 

tqpl™[n,  X  ]  =  (tcpsourceW-tqpdetW  '  ^  ^  . 

“  '  Gollimatorthickness  ^ 

/  bottcmsouroern]  -  bottatdetrn]  \  ^  ^  ^ 

bottanliJie[n,  x  ]  =  - )  x  (x-detfaoe)  +  bottcitidet[n]  ,  {n,  1,  3)1 

“  \  coUimatorthickness  '  ■' 


Plot  field  of  view  and  centerlines 


Plot[{lineA[l,x],lineB[l,x],centerline[l,x],lineA[2,x],lineB[2,x],centerline[2,x],lineA[3,x 
],lineB[3  ,x]  ,centerbne[3  ,x] } ,  {x,0,focallength}  ,PlotStyle  {  {RGBColor[  1 ,0,0] }  ,RGBColo 
r[0,0, 1  ],  {RGBColor[0, 1,0]}}  ,PlotRange  {-20,70}  ] 

(*PlotLegend  {"line  Al","line  Bl","centerlinel","line  A2","line  B2","centerline2","line 
A3", "line  B3","centerline3"},LegendPosition->{-1.5,-.l},LegendSize  {.4, .8}*) 
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200'  400  600  800  1000 


Plot[  {lineA[  1  ,x]  ,lineB  [  1  ,x],centerline[  1  ,x],topline[  1  ,x]  ,bottomline[  1  ,x]  ,lineA[2,x]  ,lineB  [ 
2,x]  ,centerline[2,x],topline[2,x]  ,bottomline[2,x]  ,lineA[3  ,x]  ,lineB  [3  ,x],centerline[3  ,x]  ,topl 
ine[3,x],bottomline[3,x]},{x,detface,sourceface},PlotStyle  {{RGBColor[l,0,0]},RGBC 
olor[0,0, 1  ],  {RGBColor[0, 1,0]},  {RGBColor[0,0,0] } ,  {RGBColor[0,0,0] } }  ,PlotLabel  - 
>"Inside  the  collimator"] 

Inside  the  collimator 


Plot  holes  on  the  collimator  front  and  back  faces 

ColFrontHoles  1  =T  able[Graphics[Rectangle[  {lineB[n,sourceface]  ,lineB  [fsourceface] } ,  {li 
neA[n,sourceface],lineA[t,sourceface] }  ]],  {n,  1 ,3 } ,  {t,  1 ,3 }  ] ; 

ColBackHoles  1  =Table[Graphics[Rectangle[  { lincB  [n,detface]  ,lineB  [t,detface] } ,  {lineA[n, 
detface],lineA[t,detface] }  ]],  {n,  1 ,3 } ,  {t,  1 ,3 }  ] ; 
ColFrontHoles2=Table[Graphics[Rectangle[{- 
hneA[n,sourceface],lineB[t,sourceface]},{- 
lincB  [n,sourceface],lineA[t,sourceface] }  ]] ,  {n,  1 ,3  } ,  {t,  1 ,3  }  ]; 
ColBackHoles2=Table[Graphics[Rectangle[{-lineA[n,detface],lineB[t,detface]},{- 
lineB[n,detface]  ,lineA[t,detface] }  ]],  {n,  1 ,3 } ,  {t,  1 ,3 }  ]; 
ColFrontHoles3=Table[Graphics[Rectangle[{-lineA[n,sourceface],- 
lineA[t,sourceface] } ,  {-lineB[n,sourceface],-lineB[t,sourceface] }  ]],  {n,  1 ,3  } ,  {t,  1 ,3 }  ]; 
ColBackHoles3=Table[Graphics[Rectangle[  {-lineA[n,detface],-lineA[t,detface] } ,  {- 
lineB[n,detface],-lineB  [t,detface] }  ]],  {n,  1 ,3 } ,  {t,  1 ,3 }  ]; 
ColFrontHoles4=Table[Graphics[Rectangle[{lineB[n,sourceface],- 
lineA[t,sourceface] } ,  {lineA[n,sourceface]  ,-lineB  [t,sourceface] }  ]],  {n,  1 ,3  } ,  {t,  1 ,3  }  ] ; 
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ColBackHoles4=Table[Graphics[Rectangle[{lineB[n,detface],- 
lineA[t,detface] } ,  {lineA[n,detface],-lineB[t,detface] }  ]] ,  {n,  1 ,3 } ,  {t,  1 ,3 }  ] ; 
Colors2=Graphics[Hue[.77]]; 

BackFace=Table[Graphics[Line[{  {-15,-1 5},  {-15,15},  {15,1 5},  {15,-1 5},  {-15,-1 5}}]]]; 
FrontFace=Table[Graphics[Line[{  {-25, -25},  {-25,25},  {25,25},  {25, -25},  {-25,-25}}]]]; 
Show[ColBackHolesl,ColBackHoles2,ColBackHoles3,ColBackHoles4,FrontFace,Colors 
2,ColFrontHolesl,ColFrontHoles2,  ColFrontHoles3,  ColFrontHoles4,AspectRatio- 
>Automatic,Axes  True] 


Plot  fields  of  view  on  the  deteetor  front  and  baek  faees 

FrontHoles  1 =T  able[Graphios  [Reetangle[  { lineA[n,orystaloenter-l-4 . 5  ]  ,lineA[t,orystaloenter 
+4.5]},{lineB[n,orystaloenter-l-4.5],lineB[t,orystaloenter-l-4.5]}]],{n,l,3},{t,l,3}]; 
BaokHolesl=Table[Graphios[Reotangle[{lineA[n,orystaloenter-4.5],lineA[t,orystaloenter- 
4 . 5  ] } ,  { lineB  [n,  ery  stale  enter-4 . 5  ] ,  lineB  [t,  erys  tale  enter-4 .5]}]],{n,l,3},{t,l,3}]; 
FrontHoles2=T  able[Graphies  [Reetangle[  { - 
lineB  [n,erystaleenter-l-4 .5  ]  ,lineA[t,erystaleenter-l-4 .5]},{- 
lineA[n,erystaleenter-l-4.5],lineB[t,erystaleenter-l-4.5] }  ]],  {n,  1 ,2} ,  {t,  1 ,3  }  ]; 
BaekHoles2=Table[Graphies[Reetangle[{-lineB[n,erystaleenter- 

4.5] ,lineA[t,erystaleenter-4.5]},{-lineA[n,erystaleenter-4.5],lineB[t,erystaleenter- 

4.5] }]],{n,l,2},{t,l,3}]; 

FrontHoles3=Table[Graphies[Reetangle[{-lineB[n,erystaleenter-l-4.5],- 
lineA[t,erystaleenter-l-4.5] } ,  {-lineA[n,erystalcenter-l-4.5],- 
lineB[t,erystaleenter-l-4.5] }  ]] ,  {n,  1 ,3 } ,  {t,  1 ,2  }  ] ; 
BaekHoles3=Table[Graphies[Reetangle[{-lineB[n,erystaleenter-4.5],- 
lineA[t,erystaleenter-4.5]},{-lineA[n,erystaleenter-4.5],-lineB[t,erystalcenter- 

4.5] }]],{n,l,3},{t,l,2}]; 

FrontHoles4=Table[Graphies  [Reetangle[  { lineA[n,erystaleenter-l-4 .5  ]  ,- 
lineB  [t,erystaleenter-l-4 . 5  ] } ,  { lineB  [n,erystaleenter-l-4 . 5  ]  ,- 
line  A[t,  cry  stale  enter-l-4 .5]}]],{n,l,2},{t,l,2}]; 
BackHoles4=Table[Graphics[Rectangle[{lineA[n,crystalcenter-4.5],- 
lineB[t,crystalcenter-4.5] } ,  {lineB[n,crystalcenter-4.5],-lineA[t,crystalcenter- 

4.5] }]],{n,l,2},{t,l,2}]; 
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Colorsl=Graphics[Hue[.77]]; 

BackStrips=Table[Graphics[Line[{{n,-12.5},{n,12.5}}]],{n,-12.5,12.5,5}]; 
FrontStrips=Table[Graphics[Line[{{-12.5,n},{12.5,n}}],],{n,-12.5,12.5,5}]; 
Show[BackHoles  1  ,BackHoles2,BackHoles3  ,BackHoles4,BackStrips,FrontStrips, Colors  1 , 
FrontHolesl,FrontHoles2,FrontHoles3,FrontHoles4,AspectRatio-> Automatic, Axes  True] 


Plot  fields  of  view  on  the  source 

SourceFIolesl=Table[Graphics[Rectangle[{lineB[n,focallength],lineB[t,focallength]},{lin 
eA[n,focallength],lineA[t,focallength] }  ]] ,  {n,  1 ,3 } ,  {t,  1 ,3 }  ] ; 

SourceFIoles2=Table[Graphics[Rectangle[{-lineA[n,focallength],lineB[t,focallength]},{- 
lineB[n,focallength],lineA[t,focallength] }  ]] ,  {n,  1 ,3  } ,  {t,  1 ,3  }  ] ; 
SourceFIoles3=Table[Graphics[Rectangle[{-lineA[n,focallength],- 
lineA[t,focallength]},{-lineB[n,focallength],-lineB[t,focallength]}]],{n,l,3},{t,l,3}]; 
SourceFIoles4=Table[Graphics[Rectangle[{lineB[n,focallength],- 
lineA[t,focallength]} ,  {lineA[n,focallength],-lineB[t,focallength] }  ]],  {n,  1 ,3  } ,  {t,  1 ,3 }  ]; 
Source=Graphics[Disk[  {0,0}  ,50]]; 

Colors3=Graphics[Hue[.77]]; 

Colors4=Graphics[Hue[.5]]; 

Show[Colors3,Source,Colors4,SourceHolesl,SourceHoles2,SourceHoles3,SourceHoles4, 
AspectRatio-> Automatic, Axes  True] 

Show[Colors4,SourceFIolesl,SourceHoles2,SourceFIoles3,SourceHoles4,Colors3, Source, 
AspectRatio-> Automatic, Axes  True] 

--  r 


Graphics 
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Graphics 
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Appendix  B,  Image  Processing  Code 

%  HPGe  Double  Sided  Strip  Detector  Imaging  code 
%  CPT  Doug  Rothenbush 
%  Version  1.0  31Dec04 

%  This  code  takes  a  text  file  of  data  from  a  HPGe  DSSD  and 
%  converts  it  into  an  imgage  of  the  source.  The  data  is  taken  from  two 
%  DGF-4C  modules.  Module  1  has  four  rear  strips,  and  module  2  has  four 
%  front  strips. 

% 

%  Code  outline 

%  1 .  Read  in  raw  data  and  user  inputs 
%  2.  Sort  the  data  into  separate  matrices  for  each  module 
%  3.  Compare  each  module  1  event  with  those  in  module  2  to  find  a 
%  two  strip  event  coincident  in  time  and  energy. 

%  4.  Determine  which  pixel  received  the  hit  and  add  1  to  its  image  value. 

%  5.  Compare  each  module  1  coincident  event  with  those  in  module  2  to  find  a 
%  three  strip  event  coincident  in  time  and  energy. 

%  6.  Determine  which  pixel  received  the  hit  (earliest  timestamp)and  add  1  to  its  image 
%  value. 

%  7.  Compare  each  module  2  coincident  event  with  those  in  module  1  to  find  a 
%  three  strip  event  coincident  in  time  and  energy. 

%  8.  Determine  which  pixel  received  the  hit  (earliest  timestamp)and  add  1  to  its  image 
%  value. 

%  9.  Display  the  image  along  with  a  surface  plot  and  a  histogram  of  hits. 

%  module  1  are  rear  strips 
%  module  2  are  front  strips 
%  1 .  Read  in  the  raw  data  and  user  inputs 
clear  %  clear  workspace 
load  'data.txf ; 

windows  =  input('How  many  energy  windows  would  you  like  to  image,  1  or  2?') 
lower_limit2  =  0; 
upper_hmit2  =  0; 
if  windows  ==  2 

lower_limit  =  input('What  is  the  lower  limit  of  your  first  imaging  window  in  keV?') 
upper_limit  =  input('What  is  the  upper  limit  of  your  first  imaging  window  in  keV?') 
lower_limit2  =  input('What  is  the  lower  limit  of  your  second  imaging  window  in 
keV?') 

upper_limit2  =  input('What  is  the  upper  limit  of  your  second  imaging  window  in 
keV?') 

range_upper_limit  =  upper_limit2; 
else 

lower_limit  =  input('What  is  the  lower  limit  of  your  imaging  window  in  keV?') 
upper  limit  =  input('What  is  the  upper  limit  of  your  imaging  window  in  keV?') 
rangeupperlimit  =  upperlimit; 
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end 

three_strip_use  =  input('Would  you  like  to  use  three  strip  events?  (1)  Yes  (2)  No 
(Default  is  yes)') 

disp('Calibration  is  of  the  form  energy  =  (channel  #  -  a)/b') 
dispCCurrently  a  =  -20.2  and  b  =  38.852') 

keep_cal  =  input('Is  this  the  calibration  you  want  to  use?  (l)yes  (2)no') 
if  keep_cal  ==  2 

a  =  input('Enter  your  new  value  for  a?') 
b  =  input('Enter  your  new  value  for  b?') 
else 
a  =  85; 
b  =  38.634; 
end 

threshold  =  b  *  lower  limit  +  a; 
threshold  max  =  b  *  upper_limit  +  a; 
threshold2  =  b  *  lower_limit2  +  a; 
threshold_max2  =  b  *  upper_limit2  +  a; 

%  initialize  image  matrices 
two_strip_image  =  zeros(4,4); 
three_strip_image  =  zeros(4,4); 
full_image=zeros(4,4); 

%  Define  time  and  energy  window  to  define  a  hit 

deltatime=10; 

deltaenergy=200; 

j=i; 

spill_size  =  data(j,l); 
modi  =  []; 
mod2  =  []; 
totalheader=[]; 
coincdatal  =  []; 
comc_data2  =  []; 
hit_hist=[]; 

%  pull  out  first  header 
header  =  data(j  :j+5 , 1 ); 
totalheader=  [totalheader  header]; 

%  jump  to  first  hit  pattern 

j=V; 

hit_patt=data(j ,  1 ); 
mod=header(2,l); 
while  mod  >  0 
%  2.  Sort  data 
while  (j  <  spill_size) 

%  find  coincident  events  for  module  1  and  put  into  coinc  datal  matrix 
if  (  hit_patt  ==  3)  I  (hit_patt  ==  6  )  %  3  =  001 1  6  =  0110 
comc_event  =  data(j;(j-l-6),l); 
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coinc_energy  =  coinc_event(5,l)  +  coinc_event(7,l); 

if  ((coinc_energy  >  threshold)  &  (coinc  energy  <  threshold_max))  |  ((coinc  energy  > 
threshold!)  &  (coinc  energy  <  threshold  max!)) 
if  mod==l 

coinc  datal  =  [coinc  datal  coinc  event]; 
elseif  mod==2 

coinc_data2  =  [coinc_data2  coinc_event]; 
end 
end 
j=j+'7; 

elseif  (  hit_patt  ==12)  %  12  =  1100 
coinc_event  =  data(j:(j+6),l); 
coinc_energy  =  coinc_event(5,l)  +  coinc_event(7,l); 

if  ((coinc_energy  >  threshold)  &  (coinc  energy  <  threshold_max))  |  ((coinc_energy  > 
threshold!)  &  (coinc_energy  <  threshold_max2)) 
if  mod==l 

coinc  datal  =  [coinc  datal  coinc  event]; 
elseif  mod==2 

coinc_data2  =  [coinc_data2  coinc_event]; 
end 
end 

j=j+7; 

else 

event  =  data(j;j+4,l); 

if  ((event(5,l)  >  threshold)  &  (event(5,l)  <  threshold_max))  |  ((event(5,l)  > 
threshold!)  &  (event(5,l)  <  threshold_max2)) 
if  mod==l 

modi  =  [modi  event]; 
elseif  mod==2 

mod!  =  [mod!  event]; 
end 
end 

j=j+5; 

end 

hit_patt  =  data(j,l); 
end 

header=data(j  :j+5 , 1 ); 
totalheader=  [totalheader  header] ; 
spill_size  =  spillsize  +  header(l,l); 
mod=header(2,l); 

%  jump  past  the  header 

j=j+6; 

hit_patt=data(j,  1); 
end 

dispCSorting  is  now  complete') 
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%  module  data  are  now  separated  and  ready  to  be  eompared 
%  3.  Determine  two  strip  events 
[modl_rows,modl_eol]  =  size(modl); 

[mod2_rows  ,  mod2_eol]  =  size(mod2); 
mmm=l; 

for  n  =  l:modl_ool 

for  m  =  minm:mod2_coi 

%  find  if  event  time  hi  words  are  equal 
if  (modl(2,n)==mod2(2,m)) 

%  find  if  fast  trigger  times  are  within  deltatime 
if  (abs(modl(4,n)  -  mod2(4,m))  <  deltatime) 

%  if  energies  are  within  deltaenergy  we  have  a  two 
%  strip  hit 

if  (abs(modl(5,n)  -  mod2(5,m))  <  deltaenergy) 
hit_energy  =  (mod2(5,m)-a)/b; 
hit  hist  =  [hit  hist  hit_energy]; 

%  4.  Determine  whieh  pixel  gets  the  hit  and  add  1  to 
%  its  image  value 

r=[]; 

f=[]; 

%  determine  whieh  rear  strip  got  hit 
if  (modl(l,n)==l) 
r=l; 

elseif  (mod  1  ( 1  ,n)==2) 
r=2; 

elseif  (mod  1  ( 1  ,n)==4) 
r=3; 

elseif  (mod  1  ( 1  ,n)==8) 
r=4; 
end 

%  determine  whieh  front  strip  got  hit 
if  (mod2(l,m)==l) 

f=i; 

elseif  (mod2(l,m)==2) 
f=2; 

elseif  (mod2(l,m)==4) 
f=3; 

elseif  (mod2(l,m)==8) 
f=4; 
end 

%  add  hit  eount  to  image  at  (f,r) 
two_strip_image(f,r)=two_strip_image(f,r)  +  1; 
minm=m; 
m=mod2_eol; 
end 
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end 

end 

end 

end 

two_strip_hits=sum(sum(two_strip_image)) 

%  end  of  two  strip  event  seetion 
%  eompute  3  strip  events 

%  5.  Compare  eaeh  module  1  eoineident  event  with  those  in  module  2  to  find  a 
%  three  strip  event  eoineident  in  time  and  energy. 
[eoine_datal_rows,eoine_datal_eol]  =  size(eoine_datal); 
forn=  Ixoinedataleol 
for  m  =  1  ;mod2_ool 

%  see  if  sum  of  eoineident  energies  are  above  the  energy  threshold 
ooino_energy  =  ooino_datal(5,n)+  ooino_datal(7,n); 

%  find  if  event  time  hi  words  are  equal 
if  (eoine_datal(2,n)==mod2(2,m)) 

%  find  if  fast  trigger  times  are  within  deltatime 
if  (abs(eoine_datal(4,n)  -  mod2(4,m))  <  deltatime) 

%  if  energies  are  within  deltaenergy  we  have  a  two 
%  strip  hit 

if  (abs(ooino_energy  -  mod2(5,m))  <  deltaenergy) 
hit_energy  =  (mod2(5,m)-a)/b; 
hit_hist  =  [hit_hist  hit  energy]; 

%  6.  Determine  whieh  pixels  reeeived  the  hit  and  add  0.5  to  their  image 
%  value. 

r=[]; 

f=[]; 

%  determine  whieh  rear  strip  got  hit 
if  (ooino_datal(l,n)==3) 

if  ooino_datal  (4,n)<ooino_data  1  (6,n) 
r=l; 
else 
r=2; 
end 

%rl=l; 

%r2=2; 

elseif  (eoine_data  1(1  ,n)==6) 

if  eoine_datal  (4,n)<eoine_datal  (6,n) 
r=2; 
else 
r=3; 
end 

%rl=2; 

%r2=3; 

else 
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if  coinc_datal  (4,n)<coinc_data  1  (6,n) 
r=3; 
else 
r=4; 
end 

%rl=3; 

%r2=4; 

end 

%  determine  whieh  front  strip  got  hit 
if  (mod2(l,m)==l) 

f=i; 

elseif  (mod2(l  ,m)==2) 
f=2; 

elseif  (mod2(l,m)==4) 
f=3; 

elseif  (mod2(l,m)==8) 
f=4; 
end 

%  add  hit  count  to  image  at  (f,r) 
three_strip_image(f,r)=three_strip_image(f,r)  +  1 ; 
%three_strip_image(f,r2)=three_strip_image(f,r2)  +  0.5; 
end 
end 
end 
end 
end 

%  7.  Compare  each  module  2  coincident  event  with  those  in  module  1  to  find  a 
%  three  strip  event  coincident  in  time  and  energy. 
[coinc_data2_rows,coinc_data2_col]  =  size(coinc_data2); 
for  n  =  1  :coinc_data2_col 
for  m  =  l;modl_col 

%  see  if  sum  of  coincident  energies  are  above  the  energy  threshold 
coinc_energy  =  coinc_data2(5,n)+  coinc_data2(7,n); 

%  find  if  event  time  hi  words  are  equal 
if  (coinc_data2(2,n)==modl(2,m)) 

%  find  if  fast  trigger  times  are  within  deltatime 
if  (abs(coinc_data2(4,n)  -  modi (4, m))  <  deltatime) 

%  if  energies  are  within  deltaenergy  we  have  a  two 
%  strip  hit 

if  (abs(coinc_energy  -  modi (5, m))  <  deltaenergy) 
hit_energy  =  (modl(5,m)-a)/b; 
hit_hist  =  [hit_hist  hit  energy]; 

%  8.  Determine  which  pixels  received  the  hit  and  add  0.5  to  their  image 
%  value. 

r=[]; 
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fi=[]; 

f2=[]; 

%  determine  whieh  front  strips  got  hit 
if  (ooino_data2(l,n)==3) 

if  ooino_data2(4,n)<ooino_data2(6,n) 

f=i; 

else 

f=2; 

end 

%fl=l; 

%f2=2; 

elseif  (eoine_data2(l,n)==6) 

if  eoine_data2(4,n)<eoine_data2(6,n) 
f=2; 
else 
f=3; 
end 

%fl=2; 

%f2=3; 

else 

if  eoine_data2(4,n)<eoine_data2(6,n) 
f=3; 
else 
f=4; 
end 

%fl=3; 

%f2=4; 

end 

%  determine  whieh  rear  strip  got  hit 
if  (modl(l,m)==l) 
r=l; 

elseif  (mod  1  ( 1  ,m)==2) 
r=2; 

elseif  (mod  1  ( 1  ,m)==4) 
r=3; 

elseif  (modl(l,m)==8) 
r=4; 
end 

%  add  hit  eount  to  image  at  (f,r) 
three_strip_image(f,r)=three_strip_image(f,r)  +  1 ; 
%three_strip_image(f2,r)=three_strip_image(f2,r)  +  0.5; 
end 
end 
end 
end 
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end 

three_strip_hits=sum(sum(three_strip_image)) 
if  three_strip_use  ==  2 

fullimage  =  two_strip_image; 
else 

full_image  =  two_strip_image  +  three_strip_image; 
end 

total_hits  =  sum(sum(full_image)) 
image_max=max(max(full_image)); 
range  =  min(hit_hist):l;max(max(hit_hist)); 

[X,Y]  =  meshgrid([l:l:4]); 

%  9.  Display  the  image  along  with  a  surface  plot  and  a  histogram  of  hits, 
if  windows  ==  2; 
subplot(2,3,[l  2  4  5]) 
surfc(X,Y,full_image) 

title( ['Surface  plot  of ',num2str(lower_limit),'-',num2str(upper_limit),'  and 

',num2str(lower_limit2),'-',num2str(upper_limit2),'  keV']) 

zlabel('Counts') 

shading  interp 

set(gca,'YDir', 'reverse') 

set(gca,'XDir', 'reverse') 

subplot(2,3,3);  imagesc(full_image,[0  max(max(full_image))]);  colormap(gray); 

title([num2str(lower_limit),'-',num2str(upper_limit),'  and  ',num2str(lower_limit2),'- 

',num2str(upper_limit2),'  keV  image']) 

xlabel('Rear  Strip  #') 

ylabel('Front  Strip  #') 

set(gca,'XDir', 'reverse') 

subplot(2,3,6),  hist(hit_hist, range) 

xlabel('Energy  [keV]') 

ylabel('Counts') 

title('Energy  Histogram') 

else 

subplot(2,3,[l  2  4  5]) 
surfc(X,Y,full_image) 

title( ['Surface  plot  of ',num2str(lower_limit),'-',num2str(upper_bmit),'  keV']) 

xlabel('Rear  Strip  #') 

ylabel('Eront  Strip  #') 

zlabel('Counts') 

shading  interp 

set(gca,'YDir', 'reverse') 

set(gca,'XDir', 'reverse') 

axis  normal 

subplot(2,3,3);  imagesc(full_image,[0  max(max(full_image))]);  colormap(gray); 
title([num2str(lower_limit),'-',num2str(upper_limit),'  keV  image']) 
xlabel('Rear  Strip  #') 
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ylabel('Front  Strip  #') 

set(gca,'XDirVreverse') 

subplot(2,3,6),  hist(hit_hist, range) 

xlabel('Energy  [keV]') 

ylabel('Counts') 

title('Energy  Histogram') 

end 
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Appendix  C,  Resolution  Measurement  Spectra 


Figure  56;  Strip  FI  activated  Sn  spectrum 


Figure  57;  Strip  F2  activated  Sn  spectrum 
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Figure  58:  Strip  F3  activated  Sn  spectrum 


Figure  59:  Strip  F4  activated  Sn  spectrum 
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Figure  60;  Strip  R1  activated  Sn  spectrum 


Figure  61;  Strip  R2  activated  Sn  spectrum 
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Figure  62;  Strip  R3  activated  Sn  spectrum 
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